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ABSTRACT 


Numerical  simulations  and  experimental  measurements  have  been  carried 
out  to  study  the  transient  behavior  of  the  P*-n”-N*  structure  under  dif- 
ferent levels  of  injection  in  both  the  reverse  and  the  forward  pulsed  direc- 
tions. A major  part  of  this  research  is  devoted  to  the  explanation  of  the 
anomalous  voltage  that  appears  across  the  diode  terminals  when  pulsed  in  the 
forward  direction  with  current  densities  in  the  range  of  10^  to  10*  arap/cm^. 

A numerical  iterative  method  of  solution  of  the  one-dimensional  basic 
two-carrier  transport  equations,  coupled  with  Maxwell's  and  Poisson's  equa- 
tions, describing  the  behavior  of  the  semiconductor  device  is  presented. 
The  method  is  of  a very  general  character,  with  none  of  the  conventional  as- 
sumptions and  restrictions  introduced,  and  freedom  is  available  in  selecting 
the  doping  profile,  generation-recombination  law,  mobility  dependencies,  and 
injection  level.  For  a specified  arbitrary  current  excitation,  the  solution 
yields  the  terminal  voltage  and  all  the  quantities  of  interest  throughout 
the  entire  device  as  functions  of  position  and  time. 

The  pure  implicit  technique  along  with  the  Newton-Raphson  iteration  are 
used  to  solve  the  equations,  while  the  explicit  technique  is  used  to  obtain 
a first  guess  to  start  the  iterative  solution. 

A numerical  method  is  arrived  at  for  the  determination  of  effective 
transient  lifetimes  in  the  case  of  excess  diffusion  due  to  very  high  current 
density.  As  the  transient  solution  approaches  steady  state,  the  transient 
lifetimes  reduce  to  their  equilibrium  values. 

The  use  of  band-to-band  Auger  recombination  as  a process  complementary 
to  the  impact  ionization  process,  along  with  the  transient  lifetimes  explain 
the  large  terminal  voltage  across  the  diode  structure  when  pulsed  with  in- 
tense current  density  in  the  forward  direction. 


CHAPTER  1 


INTRODUCTION 


Present  seal  conductor  device  Modeling  of  the  volt-aepere  characteris 


tics,  transient  response,  and  avalanche  breakdown  is  adequate  for  predict- 
ing the  response  of  junction  devices  in  the  noreal  region  of  operation.  For 


very  high  injection  levels,  such  as  diodes  operating  under  intense  current 
pulses  in  the  forward  biased  direction  with  current  density  of  10^  amp/cm^ 


these  Models  are  no  longer  adequate  and  do  not  predict  the 


anoMalous  voltage  that  appears  across  the  terainals  of  thed  device 


In  order  to  achieve  analytical  results  for  a diode  in  a dosed  fora 


solution,  a nuaber  of  assuMptions  are  usually  introduced.  These  assumptions 


liait  the  solution  to  specific  device  structures  and  operating  conditions 


Some  of  these  assuaptions  are  as  follows 


(a)  The  separation  of  the  device  into  different  regions;  either  a 


space-charge  layer  that  is  fully  depleted  of  Mobile  carriers  or 


quasi-neutral  regions 


(b)  The  neglecting  of  the  Majority  carrier  current  in  the  low  doping 


region 


(c)  The  I iaitation  of  the  doping  profiles  to  abrupt  junctions  or  to 


linearly  graded  junctions 


(d)  The  assumption  of  constant  mobilities 


(e)  The  use  of  a linear  low  level  recombination  rate 


Host  of  these  assuaptions  are  not  valid  for  high-injection  levels 


Numerical  techniques,  with  the  aid  of  high  speed  digital  computers 


represents  an  alternative  approach  to  the  problem.  This  is  accomplished  by 


solving  the  one-diaensional  two  carrier  transport  equations,  along  with 


Maxwell's  and  Poisson's  equations.  By  using  either  the  implicit  or  explicit 


- 2 - 


numerical  techniques^  one  can  obtain  the  transient  and  the  spatial  distribu- 
tions of  all  the  physical  quantities  of  interest  such  as  the  electric  field, 
the  hole  density,  the  electron  density,  the  current  densities,  etc. 

A.  De  Mari  solved  this  set  of  coupled  second  ordered  non-linear 

equations  assuming  constant  electron  and  hole  mobilities,  absence  of  a 

/ 

generation-recombination  term,  and  with  the  application  of  external  excita- 
tions of  low  to  moderately  high  current  density  values.  Scharfetter  and 
Gummel  C23  solved  a similar  set  of  equations  for  a Silicon  N*-N-P  structure 
operating  under  reversed  bias  conditions  as  a Read  diode  oscillator.  Lee, 
Lomax,  and  Haddad  C3]  developed  a computer  program  to  simulate  a Silicon 
N*-P-N^  structure  operating  as  an  avalanche  diode  oscillator.  All  of  these 
authors  used  different  numerical  techniques,  and  different  boundary  condi- 
tions to  achieve  a numerical  solution  to  the  set  of  equations. 

It  is  the  purpose  of  this  work  to  develop  a solution  to  the  set  of 
equations  that  is  free  from  any  of  the  traditional  assumptions  and  thereby 
simulate  the  transient  behavior  of  a current  driven  one-dimensional  solid 
state  structure  starting  from  the  thermal  equilibrium  condition  (or  a known 
steady-state  condition).  This  computer  simulation  is  used  to  explain  and  to 
predict  the  unexpectedly  large  transient  voltage,  which  appears  in  the  ex- 
perimental measurements  of  diodes  that  are  pulsed  with  a very  narrow  and 
large  current  pulse  in  either  the  reversed  or  the  forward  directions. 

Chapter  II  presents  the  mathematical  equations  to  be  solved  numerically 
and  the  literature  review  for  different  numerical  techniques,  along  with  the 
methods  used  in  this  research  . It  also  exposes  the  various  sources  of  er- 
rors that  could  result  in  divergence  or  instability.  Chapter  III  describes 
a numerical  method  to  obtain  an  accurate  solution  of  the  thermal  equilibrium 
pase.  The  thermal  equilibrium  solution  is  necessary  since  it  serves  as  the 
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initial  conditions  for  the  transient  problea.  As  an  exaaple  of  the  numer- 


ical calculations^  the  transient  solution  for  the  reverse  biased  case  is 
illustrated  in  Chapter  IV.  In  Chapter  the  theory  of  the  forward  biased 
case  under  different  levels  of  injection  is  discussed  in  detail  along 
with  a numerical  method  to  predict  an  effective  transit  life  time.  The 
transient  solution  is  also  presented.  Chapter  VI  describes  the  experimental 
results,  the  diode  fabrication  procedure,  and  the  method  used  to  measure  the 
physical  parameters  such  as  the  doping  profile  and  the  width  of  the  low  dop- 
ing region.  It  also  contains  a comparison  between  the  theory  and  experi- 


ment. In  the  appendices  (A  to  C)  are  gathered  the  formulation  for  the  con- 
struction of  the  matrix  coefficients,  the  algorithms  to  solve  a tridiagonal 
matrix,  and  the  tri-tridiagonal  matrix. 

I 
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MATHEMATICAL  EQUATIONS 
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It  is  the  purpose  of  this  chapter  to  present  a numerical  method  for 
the  transient  solution  of  the  one~dimensional  two  carrier  continuity  equa- 
tions coupled  with  Maxwell's  equations  and  with  Poisson's  equation.  The 
method  is  of  a very  general  character.  It  allows  for  the  freedom  of 
choosing  the  generation  recombination  laws,  the  doping  profile,  and  the  mo- 
bility dependencies.  The  transient  solution  results  in  a knowledge  of  all 
the  physical  quantities,  between  the  device  contacts,  as  functions  of  dis- 
tance and  time.  The  pure  implicit  technique  is  used  to  iteratively  solve 
for  the  electric  field,  the  electron  density,  and  the  hole  density  at 
each  instant  of  time.  The  explicit  technique  is  used  to  initiate  the 
iteration.  A literature  review  for  different  numerical  methods  is  also 
included  for  comparative  purposes. 

The  fundamental  equations  that  govern  the  one-dimensional  behavior  of  a 
semiconductor  material  are 
(i)  Maxwel I 's  Equations 

J(t)  = J^(x,t)  + Jp(x,t)  + c ^ 


3J(t)  _ Q 
ax 


(2) 


1 

i 


, 


which  shows  that  the  total  current  density  may  vary  with  time  but  is  con- 
stant with  distance.  The  third  term  of  equation  (1)  is  the  displacement 
current.  The  first  and  second  terms  are  the  electron  and  hole  currents 
respectively.  Their  current  flow  equations  are 


I 
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J„<x^t)  = e u„(x,t)  n(x,t)E{x^t)  + e l>„(x,t) 
n n n p A 


Jp(x^t)  = e Up(x^t)P(x^t)E(x,t)  - e Op(x^t)  (4) 

Th«  above  two  equations  show  that  electron  and  hole  currents  consist  of 
drift  and  diffusion  components.  The  dependence  of  mobilities  and  diffusion 

coefficients  (uCX/t)  and  D(X/t))  on  distance  and  time  follows  from  their 

\ 

dependence  on  the  electric  field  E(x,t)  as  follows 


D = — u 
P q P 


D - — u 
n q *^0 


where 


“p  = “n  = I 


v^  = 1500  E 


cm/ sec  for 


E < 3.5x10^ 


V = 1.45x10^  gO.1525  3.5x10^  < E < 1.4x10^ 

n 


v^  = 9x10“ 


cm/sec  for 


E > 1.4x10^ 


Vp  = 600  E 


cm/sec  for 


E < 3.5x10^ 


V = 4.68x10^  gO.445  for  3.5x10^  < E < 1.4x10^ 

P 


Vp  = 9x10' 


cm/ sec  for 


E > 1.4x10^ 
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(ii)  Poisson's  Equation 

■ 7 - n(x^t)  + nJ(x)  - n"(x)] 


where,  N^(x)  are  the  ionized  donor  and  acceptor  iaipurities/cn* 

(iii)  Continuity  Equations 


3n(x,t)  _ . 1 ***n^*'^^  . _ 

—if—  = - +6 ^ 


aP<x,t)  . _ , . 

at  ■ e ax 


+ 6(x,t) 


where  U(x,t)  is  the  net  hole-electron  recombination,  typically  given  as. 


U<x,t)  = 


np  - n^ 

T~7n+fi7r“+TIjf]F^ 

po  1 no  1 


If  the  life  times  are  equal,  then  = tp^,  and  the  recombination  term 
reduces  to  the  following: 


np  - n. 

U<X/t)  = X 

' T^(n+p+2n^) 


The  minority  carrier  lifetime  is  a function  of  the  doping  density  in  the 
diffused  region  and  is  therefore  shorter  than  that  of  the  bulk  region.  The 
relation  between  the  lifetime  and  the  doping  density  is  as  follows  C4J  C5D; 


where 


is  the  low  doping  region  lifetiae 


Np  is  the  doping  density 


AN1  is  a constant  typically  1.894  x 10  /ca' 


The  generation  term,  with  no  light  on  the  saaple^  is 


where 


are  the  ionization  coefficients  for  the  electrons  and  holes. 

At  this  point  it  should  be  noted  that  only  six  of  the  above  seven  equa- 
tions are  independent  since  equation  (2)  can  be  derived  from  the  remaining 


six  equations. 

With  the  application  of  a current  pulse  to  a P*-N  diode  (i.e.^  changing 
the  current  density  with  time)  the  transient  behavior  of  E(x,t),  p(x^t)^  and 
n(x,t)  becomes  of  primary  interest.  In  order  to  fully  understand  the  time 
variation  of  these  quantities,  one  must  pay  special  attention  to  the 


dielectric  relaxation  time  and  the  transit  time 


a)  The  dielectric  relaxation  time  constant 

This  quantity  cotnrol s the  variation  of  the  electric  field  with  time 
in  the  case  of  constant  electron  and  hole  densities.  An  approximate  value 
for  the  dielectric  relaxation  time  can  be  derived  if  one  assumes  that  the 
hole  and  electron  mobilities  are  independent  of  the  electric  field. 
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Co«bining  equations  (1)^  (3)^  and  (4)  for  the  displacement  current  di- 
vided by  C 

il|*4il  = lcj(t)  - J^(x)  - yx)3 


. 1 


= - CJ(t)  - e(M^n(x)  ♦ M p(x))E(x,t) 


where 


♦ e<0  - D in<l>)3 

p 9x  "n  dx 


= Y CJ(t)  - o(x)  E(x,t)  + 0(x)3 


a(x)  = e(M  n(x)  + ji„  P(x)) 
n p 


D(x)  = e(D  - D 

p 3x  “^n  ^x 


or. 


♦ -21*1  E(x^t)  = I Cj(t)  + 


D(x)] 


(8) 


If  a positive  current  step  of  magnitude  Is  applied,  the  solution  of  equa- 
tion (8)  becomes 


j - 

E(x,t)  - E(x,t=0")  = - ^ C1  - e ^ ] 


where 


E(..t. 0-)  . 


The  dielectric  relaxation  time  constant  Is  then 


L 
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_ e _ e ___ 

'4^  ” otx)  n<x)  + jip  f*c!x)J 

Therefore,  if  one  considers  a P*-N  diode  with  the  following  parameters: 

N.  = 10‘'®  * p(x) 

A 

Up  * 200 

Np  = lO**^  * n(x) 

u - 1500 
n 

G = 11.8  X 8.854  X 10"^^ 

e = 1.602  X 10"“” 

In  the  P side: 

= T-g-  = 3.261  X lO'^^sec. 
rp  eu^N^ 

and,  in  the  N side: 

T 5 — TT-  » «.347KlO"’^sec. 

'“„"d 

(2^)  The  transit  time  : 

Tiiis  quantity  controls  the  time  variation  of  the  minority  carriers  in  a 
semiconductor  region  of  width  (W),  with  a time  invariant  electric  field. 
The  general  definition  of  the  steady  state  transit  time  is: 
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^ 1 


■ J. 


where 

is  the  minority  carrier  charge  contained  in  y,  and 

J is  the  current  density  of  minority  carrier. 

m I 

The  above  transit  time  can  be  approximated  in  simple  cases: 

Ci)  For  a short  base  P*-N  diode  under  low  level  injection  conditions, 
the  minority  carrier  distribution  in  the  N region  can  be  approximated  as  in 
Fig.  1.  In  this  case. 


e W 


«n  " V ^Pn^0>  - Pno'  " « V, 


no 


is  the  diffusion  current  and  is  given  by: 


(p_(0)  - p„„) 

I _ _ ” no 

' 'p — ic 


n 


Therefore, 


= 


p 

n n ^no 

On  <Pn^0)  - p„^) 

p p n no 


n 

D 


,-4 


In  the  case  where  = 10  x 10  cm  and  D = 15.54 
n p 


,-8 


= 3.22  X 10  sec. 


(ii)  For  a P -N  diode,  with  an  N region  under  avalanche  conditions,  the 

base  region  is  under  very  high-level  injection  and  the  minority  carrier 

density  can  be  approximated  by  an  average  value  throughout  the  width  U . 

n 

The  minority  current  is  considered  as  a drift  current,  and  the  minority  car- 
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rier  velocity  is  taken  as  the  scattering  limited  saturation  velocity  tVg)- 
Therefore, 


Q - ep  U 
m n n 


J - ep  V 
m ^n  s 


The  above  is  the  ful I transit  time,  the  average  transit  time  is  about  one 
hal f this  val ue, 

W 

T = -5^  = 5 X lO'^sec. 

As  shown  from  the  above  calculations,  the  relaxation  time  is  smaller  than 
the  transit  time.  Therefore,  the  time  variation  of  the  minority  carriers 
follows  the  time  variation  of  the  electric  field. 

Mathematical  Equations  to  be  Solved 

From  the  previous  discussion,  it  is  obvious  that  the  electric  field  E 
and  the  carrier  densities  are  all  functions  of  distance  and  time.  These 
quantities  cannot  be  solved  in  the  N region  alone  because  the  boundary  con- 
ditions near  the  transition  region  have  an  unknown  time  dependence.  There- 
fore, one  must  try  to  solve  for  n,  p,  and  E for  the  entire  P^-N  device  as- 
suming the  ohmic  contacts  as  boundaries  on  both  sides  of  the  device. 

To  reduce  the  redundant  calculations,  the  parameters  of  the  basic  equa- 
tions are  normalized.  The  normalization  factors  used  here  are  the  same  as 


those  used  by  A.  De  Hari  (1).  The  following  definitions  have  been  Made: 

Dg  = unit  of  carrier  diffusion  constant  (cm^/sec) 

D (D  ) = diffusion  constant 

n p 

-e  = electron  charge 

K = Boltzman's  constant 

U = intrinsic  debye  length 


"i 
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= intrinsic  carrier  density  (=  1.18  x 10  for  Si) 

= thermal  voltage  ~ = 0.0259  volt) 

e = permitivity  (=  11.8  x 8.85  x 10”^^  F/cm  for  Si) 

u„(u„)  = mobility  of  carriers 
n p 

Table  1 lists  the  symbols,  units,  and  their  values  for  silicon  at  room  tem- 
perature. The  basic  equations  (1)  to  (7)  can  be  normalized  and  rewritten  as 
fol lows: 

Normal ized  Basic  Equations 


J_  = -v„n  - 
n n 


''n 

T 


an 

ax 


(9) 


J_  = -v_p  + 
p p 


1£ 

ax 


(10) 


||=p-n  + N^-N^ 


(11) 


an 

IT 


«KI " * ‘Mp  - - r 


(12) 
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3P  — „u  n^RW  P — ^ + P 

Tt"®''n  T^(2+n+p)  3x 


J(t)  = Jp  * •'p  " TT 

where  all  the  variables  are  now  the  normalized  variables. 

Substituting  equations  (9)  and  (10)  into  equation  (11)^  we  get  equation 
(15).  Also,  the  substitution  of  equation  (9)  into  equation  (12),  and  of 
equation  (10)  into  equation  (13)  yields  equations  (16)  and  (17). 

H ■ W n - l''nl  - “n  " - V » - •>“> 


V 2 V dV 

i"  = ^ + c - ^ ^3  If:  I7 

3 X c 


Where 


* f * * ‘"V’'  ■ -fjl-w-g) 


ap  '^P  3^P  X r V r 1 111  i£ 

Tt  ' E"  ‘•■''p  * E dE  •*  ax'*  ax 


GV  i 


-2.555x10=’ 


s 1.436  X 10^  e 
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-4.76x10^ 

B = 8.51  X 10^  e ^ 

Vn  = 37.3  E for  E < 5.11x10^ 

= 7.351  X 10^  e*1525  5.11x10^  < E < 2.2x10^  (20b) 

Vp  = 3.4  X 10^  E > 2.2x10^  (20c) 

Vp  = 12.3  E for  E < 5.11x10^  (21a) 

Vp  = 4.16  X 10^  E*^^^  for  5.11x10^  < E < 2.2x10^  (21b) 

Vp  = 3.4  x 10^  for  E > 2.2x10^  (21c) 

EQuations  (15)/  (16)/  and  (17)  are  the  three  second  ordered  coupled  non- 
linear differential  equations  that  have  to  be  solved  numerically  to  deter- 
mine E(X/t)/  n(X/t)/  and  p(X/t)  at  each  point  in  the  structure  at  each  in- 
stant of  time.  Once  this  is  accomplished/  then  the  electric  field  E(X/t) 
can  be  integrated  at  each  instant  of  time  to  obtain  the  terminal  voltage 
across  the  device  at  each  instant. 

The  boundary  conditions  are  obtained  from  the  assumption  that  perfect 
ohmic  contacts  are  made  to  the  P*  region  and  the  N region.  One  contact  is 
at  X * 0/  the  other  at  x = L.  Therefore/  the  normalized  boundary  conditions 


(19) 

(20a) 


A 
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for  the  particle  densities  are  as  follows: 


n(0,t)  = 1/N^  (22a) 

I 

p(0,t)  = (22b) 

n(L,t)  = Njj  (22c) 

i 

p(L,t)  = 1/Np  (22d) 


The  electric  field  is  chosen  at  the  boundaries  so  that  the  total  current 
equals  the  external  applied  current. 

The  initial  conditions  are  assumed  to  be  the  thermal  equilibrium  dis- 
tributions for  E(x,0)^  n(x^O),  and  p(x,0). 

Type  of  equations 

Before  discussing  the  numerical  methods  used  to  solve  the  three  coupled 
non-linear  differential  equations^  we  first  discuss  the  type  of  equations. 
The  general  form  of  a second  order  linear  partial  differential  equation  in 
two  independent  variables  x and  t may  be  written  as  follows: 


A A + 2B  ♦ c a = Ft  t> 


»x 


3t» 


at'  at' 


The  type  of  equation  is  determined  by  the  quantity: 


a = B‘  - AC 


Where  A is  the  discriminant  of  the  partial  differential  equation  (P.D.E.) 


(a)  If  A > 0,  the  equation  is  said  to  be  hyperbolic 

(b)  If  A = 0^  the  equation  is  said  to  be  parabolic 

(c)  If  A < 0^  the  equation  is  said  to  be  elliptic 

It  is  clear  that  equations  (16)  and  (17)  are  parabolic  non-linear  equations 
since  C = B 0. 

Therefore,  the  parabolic  equation  can  be  written  in  the  following  form; 


11  c 3Z 
at  - ^3^2'  37' 


Z,  X,  t) 


(23) 


To  solve  this  equation  one  needs  to  know  its  open  boundary  conditions,  i.e., 
the  value  of  Z(x,t)  has  to  be  known  on  three  sides  of  a rectangle  described 
as  fol lows: 

X = 0,  and  t ^ 0 
X = L,  and  t ^ 0 
t = 0,  and  0 < X < L 


The  obtain  a numerical  solution,  the  x and  t axes  have  to  be  divided 
into  uniform  or  non-uniform  spatial  and  time  steps  (ax  and  At).  Then  the 
spatial  and  time  derivatives  are  approximated  by  one  of  several  differential 
equation  schemes  at  each  space  and  time  point.  This  procedure  reduces  the 
partial  differential  equation  to  a set  of  algebraic  equations  to  be  solved 
simul taneousi y. 

Convergence  and  stability  are  the  two  problems  of  major  difficulty  in 
the  numerical  solution.  To  understand  what  they  mean,  the  following  terms 
are  defined: 

D = exact  solution  of  the  partial  differential  equation 
A = exact  solution  of  the  partial  difference  equation 
generated  by  the  difference  equation  scheme 
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N = nunerical  solution  achieved  from  the  partial  difference 
equation 

D - A = discretization  error 
A - N = numerical  error 

The  discretization  error  or  the  truncaton  error  is  a result  of  the  fin- 


ite sizes  of  the  time  and  spatial  steps.  The  convergence  problem  deals 
with  the  conditions  under  which  0 - A ->  0,  for  ax  and  At  ->  0.  The 
stabil ity  problem  deals  with  the  condition  under  which  A - N is  small 
throughout  the  entire  region  of  solution  (space  and  time). 

Numerical  Techniques  Used  to  Sol ve  the  Equations: 

The  normalized  set  of  equations  to  be  solved  are  as  follows: 


= J^(X,t)  + Jp(X/t)  - J(t) 


(24) 


= 6(x,t)  - R(x^t)  + 


= 6(x,t)  - R(x^t)  - 


= P(x,t)  - N(x^t)  + ND(x) 


To  solve  the  equations  numerical ly^  the  x-axis  and  the  time  axis  are 
divided  into  spatial  and  time  steps  (Ax  and  At).  The  spatial  and  the  time 
derivatives  are  then  approximated  by  a difference  scheme.  This  procedure 


reduces  the  non-linear  partial  differential  equations  to  a set  of  algebraic 
equations  to  be  solved  numerically. 

The  time-space  mesh  that  is  used  to  derive  the  difference  equations  is 
shown  in  Fig.  2.  The  one-dimensional  semiconductor  device  is  divided  into 
JW  partitions.  At  both  ends  of  each  partition^  there  are  space  charge 
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points  where  the  hole  density^  electron  density^  and  doping  density  are  de- 
fined. In  the  middle  of  the  two  space-charge  points  there  are  field  points 
I,  where  the  electric  field  intensity,  hole  current  density,  and  electron 
current  density,  and  all  the  parameters  which  are  functions  of  the  electric 
field  intensity  are  defined. 

Two  methods  could  be  used  to  transform  the  non-linear  partial  dif- 
ferential equations  into  difference  equations:  the  explicit  method  and  the 
implicit  method. 

In  the  explicit  method,  if  all  the  variables  are  known  at  any  time  lev- 
el t|^,  then  the  variables  can  be  calculated  directly  (that  is  explicitly)  at 
the  time  level  t|^^^  by  evaluating  the  spatial  derivatives  at  time  level  t|^. 
The  implicit  method  is  a more  complicated  procedure.  It  consists  of 
representing  the  spatial  derivatives  by  a finite-difference  form  evaluated 
at  the  advanced  time  level  t|^^^  instead  of  t|^  as  in  the  explicit  method. 
This  results  in  a system  of  simultaneous  nonlinear  equations  in  the  unknown 
variables  at  time  level  The  pure  implicit  technique  is  used  here, 
with  the  aid  of  the  Newton-Raphson  iteration  method  to  solve  the  equations, 
while  the  explicit  technique  is  used  to  calculate  a first  guess  for  the 
unknown  variables  in  order  to  start  the  Newton-Raphson  iteration.  The  basic 
reasons  for  this  choice  will  be  discussed  in  a later  section  of  the  report. 


The  pure  impi icit  technique: 

The  time-space  mesh  of  Fig.  2 is  used  to  convert  the  partial  differen- 
tials into  algebraic  equations  in  a difference  equation  form.  It  is  assumed 
that  the  three  variables  of  interest,  (the  electric  field,  the  hole  density, 
and  the  electron  density)  are  known  at  time  level  t|^  and  at  all  the  previous 
time  levels.  In  the  rest  of  this  chapter,  the  time  index  is  assumed  to  be 
K 1,  unless  otherwise  specified.  For  example,  E(I)  means  the  electric 


Hlik 
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field  at  point  (1/  K + 1)/  while  E^(I)  means  the  electric  field  at  point 
(I/K).  Keeping  this  notation  in  mind^  the  different  quantities  in  equations 
(24),  (25),  and  (26)  are  evaluated  at  points  (J,  K + 1)  and  (I,  K + 1)  as 
follows: 


3E(I) 

_ Ed) 

- e'^(I) 

3t 

At 

3P(J) 

P(J) 

- P'^(J) 

3t 

At 

3N(J) 

_ N(J) 

- n'^(j) 

3t 

At 

(28) 


The  drift  components  of  the  hole  and  electron  current  densities  at 
point  (I)  depend  on  whether  the  electric  field  is  positive  or  negative  as 
explained  in  Fig.  3. 

For  negative  electric  field,  the  drift  current  components  depend  on 
N(J)  and  P(J2),  while  in  the  case  of  a positive  electric  field  they  depend 
on  P(J)  and  N(J2).  Using  a central  difference  for  the  diffusion  current 
density  terms,  the  difference  equations  for  the  current  densities  become: 

For  Ed)  > 0 


JN(I)  = -VN(I)*N(J2)  - CN(J2)  - N(J)3 


JP(I)  = -VP(I)*P(J)  + 

For  E(I)  < 0 


|VP(I)|  1 


|E(I) 


-i-  CP(J2)  - P(J)] 
AX 


(29) 


JN(I)  = -VN(I)*N(J)  - CN(J2)  - N(J)] 


JP(I)  = -VP(I)*P(J2)  ^ CP(J2)  - P(J)D  (30) 

The  generation  term  at  the  Jth  point  is  approximated  by  its  average  value 
between  the  Ith  and  the  (I  - 1)th  point  (I  and  ID  as  shown  in  Figure  2. 

G(J)  = j (G(ID  + G(I)) 

= 2 Ca(ID*|JN(lD|  + B(ID»|JP<I1)n 


yhile 


+ Co(I)*|JN(I)|  + B(I)*|JP(I)|3 


R(J) 


N(J)*P(J)-1 
TgC2+N(J)  + P(J)3 


(3D 


(32) 


— - JPCIDD  (33) 

^•*a!**-  = t CJN(I)  - JNdDD  (34) 

dX  AX 

It  is  clear  that  the  quantities  in  equations  (28)  through  (34)  depend 
directly  or  indirectly  on  the  electric  fields  the  electron,  and  the  hole 
densities,  all  of  which  are  defined  at  points  Jl,  II,  J,  I,  J2,  and  12. 
Therefore,  the  vector  of  unknowns  in  equations  (24),  (25),  and  (26)  at 
points  I and  J is  defined  as  follows: 

CX]  = CE(ID,N(JD,P(JD,E(I),N(J),P(J),E(I2),N(J2),P(J2)D^  (35) 

By  writing  equations  (24)  to  (26)  in  a difference  form,  and  by  moving  all 
terms  in  the  three  discrete  equations  to  the  right  hand  side,  and  by  using 
equations  (28)  through  (30)  to  show  the  dependency  of  each  term  on  the  vec- 


tor  we  obtain  the  following  three  discrete  functions: 


F1I  = liiLziiiLl  + JN(I)  CE(I)^N(J)^N(J2)] 

aZ 

+ JP(I)  CE(I),P(J),P(J2)J  - J(t)  = 0 (36) 

F2J  = + 1 6(I1)CE(I1)^N<J1)^P(J1)^N(J)^P(J)D 

At  c 

+ y 6(I)CE(I),N(J),P(J)^N(J2),P(J2)]  - R(J)CN(J)^P(J)] 

+ <JP<I>CE(I>rP<J>/P^J2)3  - JP(I1)CE(I1),P(J1),P(J)]>  (37) 

AX 

F3J  = n'^(J)  - N(J)  ^ 1 g(I1)CE(I1)^N(J1)^P(J1),N(J),P(J)] 

At  c 

+ ^ 6(I)CE(I)^N(J)/P(J)^N(J2),P(J2)3  - R(J)CN(J),P(J)D 
- 1 {jN(I)[E(I)^N(J)^N(J2)]  - JN(I1)CE(I1),N(J1),N(J)]>  (38) 

AX 

The  one-dimensional  semiconductor  device  contains  a total  number  of  "M" 
spage-charge  points,  and  "M-1”  field  points,  including  the  boundaries  J = 1 
and  J = n. 

At  time  level  K + 1,  equation  (36)  will  be  written  once  for  each  point 
2£IjCH-1,  while  equations  (37)  and  (38)  will  be  written  for  each  point 
2£J£H-1,  resulting  in  a system  of  3(N-1)  simultaneous  equations  in  the 
3(H-2)  unknowns. 
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E(I)^N(J)/P(J),  where  2 < I,  J < M-1 . 


The  boundary  conditions  for  the  particle  densities  are  obtained  from 
the  assumption  that  perfect  ohmic  contacts  are  made  to  the  device  termi- 
nals. One  contact  is  at  x = 0/  the  other  is  at  x = L.  Therefore,  the  nor- 
malized particle  densities  at  the  boundary  are: 

N(0,t)  = N(1)  = 1/NA’ 

P(0,t)  = P(1)  = na" 

(39) 

N(L,t)  = N(M)  = Njj 

P(L,t)  = P(M>  = 1/ND* 

The  electric  field  is  obtained  iteratively  through  the  solution  of 
equations  (36)  to  insure  a constant  total  current  at  the  boundaries,  and 
throughout  the  whole  device. 

Newton-Raphson  Iteration  Method 

To  solve  the  system  of  nonlinear  3(M-1)  simultaneous  equations,  the 
Newton-Raphson  iteration  method  was  used.  In  order  to  have  a better  under- 
standing of  this  procedure,  Newton's  Method  for  solving  a single  non-linear 
equation  is  discussed  first. 

This  method  is  illustrated  graphically  in  Fig.  4.  If  this  figure,  if 
F(x)  = 0,  then  x = a is  the  required  solution.  If  an  initial  approximation 
x^  to  the  root  a is  provided,  a sequence  he  defined  by 

the  recursion  relation 


3F(x)  I 
~Si  'x=x^ 


F(x.) 

’‘i"*i+1 


or 


*1+3  1+2  1+1 

Figure  4,  Newton's  Method  to  solve  F(x)  ■ 0 


(40) 


F'(x^)  • - x^]  = “F(x^) 


We  draw  a tangent  to  the  curve  y = F(x)  at  the  point  (x^,F(x^)).  This 


tangent  meets  the  x-axis  at  the  point  (x.  ,^0).  If,  then,  the  curve  crosses 

L'' • I :•  , ^ ' * ■ . ^ ~ * 

the  x-axis  at  a point  (a,0)  sufficiently  near  (x.,F(x.)),  it  may  easily  be 
seen  that  the  number  x^^^  is  nearer  o than  x^.  A criterion  for  terminating 


the  iterative  procedure  of  (40)  is: 


1 ~ ^i 
’‘i  + 1 


Here  c is  a small  positive  number.  For  e = 10"",  the  final  value  of  x^^^ 
should  be  accurate  to  approximately  N significant  figures.  It  should  also 
be  noted  that  Newton's  method  may  fail  to  converge  to  a root,  if  given  a bad 
starting  value. 

The  concept  behind  the  Newton-Raphson  Hethod  is  identical  except  that 
it  applies  to  solving  a set  of  simultaneous  nonlinear  equations.  Therefore, 
equation  (40)  will  take  the  form  of  a matrix  equation  as  follows: 


(x''>  - x'-l^  = -F(x''-'^ 


where  x^^  is  the  old  approximation  of  the  3(M-2)  arrays  of  unknowns, 
while  x^^^  is  the  better  approximation. 

F(x^^”^^)  is  an  array  of  3(M-2)  functions  similar  to  equations  (36), 


(37),  and  (38)  evaluated  according  to  the  old  approximation  x 


(i-1) 


r 


p 


i 


F^^(X)  = 


3X, 
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where  2 ^ L £ H-1 
and  2 < K < H-1 


(42) 


With  the  above  notation  in  mind,  equation  (41)  is  applied  to  equations 
(36),  (37),  and  (3E)  to  evaluate  the  corrections  in  the  array  CXD  of  equa- 
tion (35).  Therefore,  the  discretized  equations  at  points  I and  J yield  the 
following  matrix  equations: 


1 -2  “3l  ‘’l  ^*2  ‘^3|'l  *^2  *^3 


a. 

a. 


a,i 

*6 1 ^ 


b. 

b. 


’6^4 


AE(I1)' 

r *1 

AN(J1) 

AP(J1) 

-FIKX^-^) 

AE(I) 

AN(J) 

2 

-F2J  (X^"'') 

AP(J) 

AE(12) 

AN(J2) 

-F3J(X^"‘') 

AP(J2) 

- 

(43) 


®7  ®8  ®9 1 ^7  ‘*9  I ^7  *^9] 


Equations  (36),  (37),  and  (38)  show  the  dependency  of  F1I,  F2J,  and  F3J 
on  the  elements  of  the  X array.  Therefore,  the  Matrix  elements  in  equation 
(43)  are  defined  as  follows: 


_ 3F1I 
®1  ■ 3E(I1) 


= 0 


- _ 3F1I  _ n 
®2  ° 


- 3F1I  _ 


®3  ' 3P(J1) 


= 0 


a,  = 


3F2J  _ 1 36(11) 
3E(I1)  ■ 7 3E^11) 


3F2J  _1  36(11) 
3N(J1T  3NCJTT 


1 3JP(I1) 

AX  3E(I1) 


(44a) 


(44b) 


(440 


(44d) 


(44e) 


. aF2J  _ 1 3G(I1)  _ 1 ajp(il) 
6 ■ 3P(J1)  ■ 7 3P(J1)  " Ax  3P(J1) 


<44f) 


_ 3F3J  _1  36(11) 
7 " 3E(I1)  1 3E(I1) 


1 3JN(I1} 

Ax  3E(I1) 


1 3JN(11) 

Ax  3N(J1} 


_ 3F3J 
8 ■ 3N(J1) 


1 36(11) 

2 3N(J1) 


(44h) 


, _ 3F3J  _ 1 36(11) 

'9  ■ wmr  ■ 7 3P(ji) 


(44i) 


_ 3F1I  _ -1  . 3JN(I)  . 3JP(1) 
1 " 3E(I)  ■ At  3E(I)  3E(I) 


(45a) 


_ 3F1I  . 3JN(I) 
2 ■ TRTjT  ■ 3N(J) 


(45b) 


. 3 FI  I 
'3  ■ 3P(J) 


3JP(1) 

3P(J) 


(450 


_ 3F2J  _ 1 36(1) 
4 ■ TlTn  ■ 7 3E(I) 


1 3JP(1) 

AX  3E(I) 


(45d) 


_ 3F2J  . 1 36(11)  , 1 36(1)  _ 3R(J) 
5 ■ TnUT  “ 7 3M(j)  * 7 WUT  WUT 


(45e) 


, - 

6 ■ 3P(J) 


1 36(1) 
7 3P(J) 


zl  ♦ 1 36(11) 
At  7 3P(J) 


3R(J) 

3P(J) 


(45f) 


_ 3F3J  _ 1 36(1)  _ 1 3JN(1) 
7 ■ 36(1)  ■ 7 3E(l)  " Ax  3E(l) 


1 36(1) 
7 3N(J) 


_ 3F3J  _ -1  . 1 36(11) 
8 ■ 3N(J)  "At  7 3N(J) 


3R(J) 

3N(J) 


(45h) 
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. _ »F3J  . n 

*=7  ■ Timr  ■ ° 


aF3J  _1  36(1)  1 3JN(I) 

*^8  ■ 3N(J2)  T aN(J2)  “ Ax  aN(J2) 


3F3J  1 36(1) 

*^9  ■ aP(J2)  ■ 7 3P<J2) 

A detailed  calculation  of  these  elements  is  found  in 
Poisson's  Cequation  (27)3  is  used  to  eliminate  the  AE(I1) 
a^^  a^/  and 

The  discretized  form  of  Poissons' s equation  is: 


L. 


(45i) 

(46a) 

(46b) 

(46c) 

(46d) 

(46e) 

(46f) 

(46g) 

(46h) 

(46i) 

Appendix  A. 
coefficients. 
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(47) 


If  for  any  reason^  the  values  of  N(J)  and  P(J)  were  changed  by  AN(J)  and 
AP(J)^  the  changes  of  AE(I)  and  AE(I1)  will  be  related  by  equation  (48) 


AE(I1)  = AE(I)  - Ax  CAP(J)  - AN(J)3 


(48) 


Therefore^  the  coefficient  AE(I1)  will  be  distributed  to  the  coeffi- 
cients of  AE(I),  AP(J)^  and  AN(J)  as  follows: 


a^AEdD  = a^AE(I)  - (a^Ax)  CAP(J)  - AN(J)] 


a^  AE(I1)  a^AE(I)  - (a^Ax)  CAP(J)  - AN(J)3 


This  will  change  the  elements^  a^,  a^,  b^,  b^,  b^,  b^,  bg,  b^  as  fol- 


lows: 


a^  - 0 


a^  ♦ 0 


bj  bj  + (a^Ax) 


b^  -»  b^  a^ 


(49a) 


(49b) 


(49c) 


(49d) 


(49e) 


(49f) 


bg  ♦ bg  + (a^Ax) 


(49g) 


b^  ♦ b^  - (a^Ax) 


(49h) 


Rewriting  equation  (43)  using  equations  (44)^  (45),  (46),  and  (49)  results 


in  equation  (50). 


0 0 0 I b^  b2  bjio  C2  Cg 

0 aj  ag|  b^  bj  bg  | 0 Cj  Cg 

0 ag  a^i  b^  bg  b^  | 0 Cg  c. 


■AE(I1)‘ 

AN(J1) 

AP(J1) 

AE(I) 

AN(J) 

AP(J) 

AE(I2) 

AN(J2) 

AP(J2) 


-FIKX’"^) 


-F2J(X^"^) 


-F3J(x'"S 


By  defining  the  subiaatrix  CA.],  CB^3,  and  CC^3,  as  follows. 


0 0 0 


bf  b2  bj 


0 C2  Cj 


CA.]  = 0 aj  ag  , CB^3  = b^  b^  bg  , CC.D  = 0 


*=5  *=6 


° *8  *9! 


'■b  '>91 


" U '9 


and  by  repeating  equation  (50)  for  all  the  spatial  points  2 £ I ^ M-1,  and 


I 


I n 


2 < J < H-1/  at  tine  level  of  K + we  obtain  the  following  equation 


AE(2) 

AN(2) 

AP(2) 


AE(I) 

AN(J) 

AP(J) 


AE(H-1) 

AN(M-1) 

AP(H-1) 


where  CJ]  is  a Tri-tridiagonal  Matrix  with  a dimension  of  3(N-2) . Also  CAp] 


and  were  eliminated  due  to  a knowledge  of  boundary  conditions 


Equation  (52)  is  solved  at  each  Newton-Raphson  iteration.  The  algo 


rithm  of  solving  this  equation  is  included  in  Appendix  B C93 


The  solution  of  equation  (52)  gives  CAx'’  ^3/  which  is  the  difference 


between  the  two  iterations  i and  i-1.  When  this  difference  reaches  a very 


small  number,  we  consider  that  the  solution  for  the  instant  of  time  t 


satisfactory.  Using  this  solution^  we  proceed  to  solve  for  the  next  instant 
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[ To  start  the  Iteration  procedure  at  any  instant  of  tiae,  a first  esti- 

aate  for  the  array  of  unknowns  is  calculated  using  the  explicit  technique. 

The  expl icit  technique. 

The  space-tine  nesh  used  for  the  explicit  technique  is  shown  in  Fig.  5. 
In  this  nethod^  the  right  hand  side  of  equations  (24)  to  (26)  are  evaluated 

at  the  time  level  t|^.  The  definitions  of  JN,  JP,  G,  and  R are  the  same  as 

in  the  implicit  method.  This  allows  us  to  directly  calculate  the  different 
variables  at  time  level  from  a knowledge  of  their  values  at  the  pre- 

vious time  level  t|^.  Therefore^  Maxwell's  equation  and  the  two  continuity 
equations  for  points  1 and  can  be  written  as: 

£*^*^(1)  * e'^(I)  + AtCJN'^(I)  + jp'^(l)  - j'^(t)D  (54) 

p’^^^J)  = p’^(J)  ♦ AtC6*^(J)  - R*^(J)  (jp'^(l)  - Jp'^(I1))T  (55) 

n'^^^(J)  = n'^(J)  + AtC6'^(J)  - R*^(J)  +1^  (Jn'^(I)  - Jn'^(I1))]  (56) 

It  is  clear  that  the  value  of  E(I)  at  time  level  tj^^^  depends  only  on 
the  values  of  N(J)^  P(J)^  E(I),  N(J2),  and  P(J2)  at*the  level  t^^,  while  the 

values  of  N(J)  and  P(J)  at  the  level  tj^^^  depends  on  N(J1)^  P(J1),  E(i'1), 

N(J)^  P(J)^  E(I),  N(J2),  and  P(J2)  at  time  level  t^^. 

Referring  to  Fig.  t,  only  those  variables  within  the  pyramid-shaped 
area  A can  have  any  influence  on  the  value  of  E(I)^  P(J),  N(J)  at  time  level 
tK+V  known  that  the  solution  of  the  partial  differential 

equations  depends  also  on  the  variable  values  in  B for  times  earlier  than 

tR+v  Therefore^  this  method  is  used  here  only  to  calculate  a first  guess 
of  the  variables  to  start  the  Newton-Raphson  iteration  at  time  level  t^  . . 
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To  start  the  transient  solution,  an  accurate  initial  condition  is  also 
required  as  a boundary  condition,  which  describes  the  spatial  distribution 
of  the  electric  field,  the  hole  density,  and  the  electron  density  at  time  t 
= 0.  The  initial  conditions  used  here  are  the  thermal  equilibrium  condi- 
tions, Chapter  III  explains  in  detail  the  method  used  to  obtain  the  thermal 
equilibrium  conditions  along  with  some  computer  generated  results. 

Literature  review  for  different  numerical  techniques 

A.  De  Mari  C13  solved  a similar  set  of  coupled  non-linear  differential 
equations  using  the  following  assumptions: 

(i)  constant  hole  and  electron  mobilities,  i.e.,  the  region  where 
the  hole  and  electron  velocities  are  linearly  dependent  on  the 
electric  field,  was  considered. 


(ii)  the  absence  of  the  recombination-generation  term. 


(iii)  external  current  excitation  in  the  form  of  a current  step  with 
a maximum  value  of  current  density  smaller  than  what  is  con- 
sidered for  this  research. 

Using  the  generalized  pure  implicit  scheme,  equation  (23)  was  represented  by 
the  following  difference  equation: 


= I^F.  = F <^tci^i^3x.,^tcH^Dx.,z(x,V> 


where 


K is  the  index  identifying  the  instant  of  time  tj^ 
i is  the  index  identifying  the  spatial  point  x 


is  the  magnitude  of  the  time  step  at  the  instant  K 


In  DeHari's  work  a non-uniform  two-dimensional  mesh  in  the  position  and 


time  coordinates,  as  shown  in  Fig.  7,  was  used.  The  assigned  boundary  and 
initial  conditions  specify  al I the  quantities  at  i = 1,  i = L,  and  K = 0. 


The  boundary  conditions  used  for  the  particle  density  were  the  same  as 


presented  in  equation  (22).  As  for  the  electric  field,  charge  neutrality 


was  assumed  at  the  ohmic  contacts.  This  boundary  condition  can  be  derived 


from  equation  (5)  to  give  the  following  equation 


The  initial  conditions  for  n,  p and  E were  either  the  thermal  equilibrium 


distributions  or  an  arbitrary  steady  state  distribution.  Therefore,  the 


three  non-linear  equations  generate  the  following  system  of  difference  equa 


tions 


(FID 


(Fill) 


Quantities  at  instant  of  time  t < t,._,,  at  any  spatial  point  0 < x.  < L are 


assumed  known.  An  iterative  procedure  is  employed  at  each  instant  of  time 


to  determine  the  unknown  quantities.  Therefore,  a third  index  j was  intro 


duced  to  identify  the  cycle  of  iteration,  while  the  correction  for  the  elec 
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tron  and  hole  densities  An  and  Ap  were  defined  as  the  difference  between 
densities  corresponding  to  two  successive  iterations. 


_ „Cj+1)  _ „(j) 
K^"i  - K"i  K"i 


(61a) 


(j+1)  = (j+1)  . (j) 

K k^i 


(61b) 


The  quantities  labeled  K in  equation  (59)  and  (60)  were  identified  with  the 
subscript  (j+1)  and  by  using  relations  (61)^  the  following  equations  were 


obtained: 


K n K K n 


= ,.(Fr)^^^  (64) 

K a K E a 

where  and  Tp  are  tridiagonal  matrices.  The  quantities  An^  Ap,  and  Fp 
are  the  vectors: 


I 3 
r 


r 
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<^>2 


(F  ).  . 
n 1-1 


^V3 


'Vi-1 


Sin9e  the  iterative  procedure  is  repeated  at  each  instant  t,  the  index  j was 
set  to  zero  each  time  consistency  was  reached  at  one  instani.  Suantities 
labeled  with  the  subscript  j = 0 are  then  the  trial  distribution  required  to 
initiate  the  iterative  procedure  at  the  next  instant  of  tine. 

Scharfetter  and  Gunnel  C2]  solved  a sinilar  set  of  equations  for  a sil- 
icon Read  diode  oscillator.  Their  equations  did  not  directly  contain  a dis- 
placenent  current  conponent.  Therefore,  after  using  the  difference  approxi- 
nation  in  Poisson's  equation  and  the  continuity  equations,  the  following  set 
of  equations  were  obtained: 


= g(N)  - CJ  (H)  - J (H-1):/Ax 


= g(N)  + - J^(H-1)]/Ax 


— ^ Cp(N)  - n(N)  + N.(N)  - N.(N)]  (67) 

Ax  c DA 

where  the  Mth  point  is  located  nidway  between  the  najor  nesh  points  N 1 
and  N.  The  current  equations  were  taken  as: 


Jp  = q .pPE  - KT  pp  ^ 


(68a) 


J„  = q M„nE  ♦ KT  tr 


(68b) 


It  was  noted  that  substituting  expressions  (68a)  and  (68b)  in  (65)  and  (66) 


leads  to  a nunerical  Instability  whenever  the  voltage  change  between  mesh 
points  exceeds  2KT/q.  Therefore/  equations  (68a)  and  (68b)  were  treated  as 
differential  equations  in  p and  n with  p^/  and  E assumed  con- 

stant between  mesh  points.  The  solution  of  these  then  relates  J and  J to 

n p 

the  other  variables  as  follows: 


Jp(H)  = E(H) 


P(M)Pp(M) 


Jj^(H)  = E(H) 


n(N+1)p^(H) 

l-exp(-E(M)Ax)) 


P(N+l)p^(H) 


n(N)p^(N) 

l-exp(E(N)Ax)) 


(69a) 


(69b) 


The  hole  and  electron  concentrations  were  presented  by  the  two  component 
vector  Y;  where/ 


0 

Y = f(Y) 


i.e./  the  time  derivative  of  Y equals  some  nonlinear  function  of  Y.  If  Yg 
is  the  initial  value  of  Y/  and  aY  is  the  deviation  from  this  value  as  time 
increase's;  and/ 


y(t)  = Yn  * ay(t) 


«y  = ftyn)  ♦ N«y  = 


ill  +ii  6y 
^ t=0 


If  M/  f/  and  SY  were  ordinary  numbers/  the  solution  of  (44)  is: 


«y(t)  = (e"^-1)  f/n 
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6y(t)  = f«t(1  + Ht/2  + + ...) 


For  small  but  finite  time  step  this  can  be  written  in  the  following  form: 


«y(t)  = 


ft 


1 - 


Nt 


or 


(1  - ^)  fiy  = ft 


The  vector  Y was  defined  as:  Y^  = 

n n 


(71) 


where  y^,  y^,  y^  are  deviations  in 


the  hole  concentration^  electron  concentration,  and  electric  field.  The 
right  hand  side  of  the  continuity  equation  is  defined  as  the  vector  f^^  as 


fol  lows;  fjj  = 
of  time  is; 


The  vector  equation  that  was  solved  for  each  increment 


(y,)  . 


(72) 


where 


3f. 


K.  = 


N 3y, 

The  matrix  elements  of  M are  defined  by  the  expanded  form  of  (72)  shown 
below: 


N-1  ^ 


+ CXR] 


N+1 


(73) 


♦ (FL)p,  yj  + (FR)|^  Xj  = (RH)^^ 
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where:  (RH)|^  = 2(fQ)|^,  the  2x2  matrices  XL^  XM,  and  XR  are  obtained  by 
taking  the  partial  derivatives  of  with  respect  to  y^  and  y2  at  mesh 
points  N-1^  N,  and  N+1  respectively^  while  the  two  component  vectors  FL  and 
FR  are  obtained  by  taking  the  partial  derivatives  of  f^^  with  respect  to  y^ 
at  mesh  points  M and  M-1. 
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The  boundary  conditions  were  Introduced  by  requiring  that  In  equa- 
tion (68b)  Is  equal  to  the  terminal  current  density  at  the  N-contact  and 
that  Jp  In  equation  (68a)  Is  equal  to  the  terminal  current  density  at  the 
P-contact.  In  addition^  the  electric  field  at  the  two  end  metallic  contacts 
Is  assumed  to  be  zero. 

Lee,  Lomax,  and  Haddad  C3]  also  solved  a similar  set  of  equations  for 
SI  avalanche-diode  oscIMtors,  using  the  Implicit  and  explicit  techniques. 
The  recombination  term  was  Ignored  In  using  the  Implicit  technique.  In  the 
explicit  method,  the  time  space  mesh  shown  In  Fig.  5 was  used  to  derive  the 
difference  equations.  The  unknown  parameters  at  the  Instant  t = t^^^^  were 
obtained  as  follows: 


j 


1 

i 

i 


„K+1  _ , r9Pi 

Pj  - P.  + At  • ^ ^ 


= Nj  + At  • _ 


,K+1  . pK+1  . riil 
El  - Ei-1  ♦ AX  L— ^ 


K+1 


For  each  time  advancement,  the  time  step  At  Is  not  a constant.  It  has  to  be 
determined  at  each  time  step  to  satisfy  certain  numerical  stability  condi- 
tions In  order  to  Insure  that  the  numerical  model  does  not  result  in  oscil- 
latory behavior  and  diverge. 


In  the  Implicit  method,  the  time  space  mesh  shown  in  Fig.  2 was  used 
The  vector  Q and  the  vector  F were  defined  by: 


The  objective  was  to  find  a difference  equation  to  approximate  the  continui 


ty  equation  evaluated  at  the  instant  of  time  halfway  between  the  Kth  time 


step  and  the  (K+1)th  time  step  or 


dQ,K+1/2 

at-* 


K+1/2 


Equation  (75)  is  identical  to  equation  (73)  of  Scharfetter  and  Gummel 


although  the  argument  used  in  deriving  equation  (75)  was  different.  The 


above  equation  is  solved  for  each  instant  of  time 


At  the  boundaries,  it  was  assumed  that  the  surface  state  density  is 


high  enough  so  that  surface  recombination  velocity  is  infinite  and  the  sem 


iconductor  is  extrinsic;  hence  the  particle  densities  at  the  boundaries  as 


sumed  their  thermal  equilibrium  values.  For  D.  C.  solutions  with  constant 


current,  the  value  of  the  electric  field  at  the  boundaries  was  chosen  itera 


tively  such  that  the  space  average  of  particle  density  over  the  active  re 


gion  of  the  device  was  equal  to  the  specified  D.  C.  current 


FP  _ 

3P/3t 

FN 

3N/3t 

As  a conclusion,  the  explicit  method  is  very  inexpensive  for  each  ad- 
vancement of  At  but  because  the  size  of  At  is  limited  by  the  stability  con- 
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ditions/  it  may  become  expensive  when  the  particle  densities  are  high^  the 
mobilities  are  high,  or  when  the  simulation  is  to  be  carried  over  a long 
time  duration.  The  implicit  method  is  complicated  and  is  more  expensive  for 
each  At,  but  it  is  more  stable  numerically.  Also,  the  time  step  At  can  usu- 
ally be  made  very  large.  It  has  been  observed  that  if  the  initial  guess  is 
not  good  the  implicit  method  sometimes  becomes  unstable  when  A t equals  in- 
finity, i.e.  for  steady  state  calculations. 
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CHAPTER  III 

THERMAL  EQUILIBRIUM  CONDITION 

This  chapter  describes  a nmierical  Method  used  to  obtain  an  accurate 
distribution  for  the  electric  field,  the  hole  density,  and  the  electron  den- 
sity at  thermal  equilibrium.  The  thermal  equilibrium  solution  is  used  as  a 
boundary  condition  for  the  transient  problem  as  discussed  previously  in 
Chapter  II. 

In  this  formulation  one  has  complete  freedom  of  choosing  the  semicon- 
ductor device  structure  between  the  two  ohmic  contacts.  Thi'-  's  accom- 
plished by  specifying  the  acceptor  and  the  donor  impurity  distribution. 
Although  the  formulation  presented  could  easily  be  modified  to  include  de- 
generate semiconductors,  Boltzmann  statistics  are  assumed. 

The  thermal  equilibrium  condition  is,  by  definition,  with  no  applied 
voltage  and  no  net  current  flow.  The  solution  is  obtained  by  accurately 
solving  Poisson's  equation: 


= I [P(x)  - N(x)  + nJ(x)  - N‘(x)] 


(76) 


ir 

L 


A semiconductor  device  consists  of  bulk  regions  and  transition  regions, 


As  an  example.  Figure  8 shows  the  electron  energy  band  diagram,  and  the  po- 

+ 


tential  variations  with  distance  for  a P -N  structure.  By  using  the  follow- 
ing definitions: 


♦j.(x)  = V(x)  = Ep(x)  - E.(x) 


and. 


(77a) 


Band  diagram  and  potential  variation  for 
structure 

(a)  electron  band  diagram 

(b)  potential  variation  with  distance 


(77b) 


(77c) 


Each  term  of  equation  (76)  is  written  as  a function  of  ^(x)  as  follows 


where  and  E-  are  the  donor  level  and  the  acceptor  level  as  shown  in  Fig 


(83a) 


(83b) 


Therefore  equations  (81)  and  (82)  can  be  rewritten  as 


Donor  and  acceptor  levels 

(a)  for  N-type 

(b)  for  P-type 


Figure  9 
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No(x)  = Np(x) 


N^(x)  = N^(x) 


[i  ♦ 2a 

[«  t bV*'*’] 


where 


a = e 


W'S  - 


b = e 


7T 


(E.  - Eg/2) 


By  substituting  equations  (78)  to  (80)^  and  (84)  and  (85)  into  equation 
(76),  equation  (86)  results: 


J^4(x) 


2 (€*VT) 


n..-*'-’  - 

L ' ’ 1 * 2 a e*'"’ 


1 + 2 a e’ 


1 

4 ♦ b 


In  order  to  solve  equation  (86)  nunerically,  the  device  is  divided  into 
spatial  steps.  The  second  order  derivative  is  approximated  by  the 
three-point  finite  difference  formula.  Therefore,  for  point  J in  Fig.  10, 
equation  (86)  reduces  to  the  following: 


>(J2)  - 4(J)  4(J)  - a(Jl) 

H(J)  “ H(J1) 


Tc"?VTy 
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1+2ae’ 


^(J)  1 

.-♦(J)j 


GC(KJ1)^«(J)/<I«CJ2)D  = 0 


= H<J)(|.(J1)  - (H(J)+H(J1))*(J)+H{J1)^(J2) 


. qH(J)H(Jl)(H(J)+H(J1)) 
2.e.VT 


r /.X  N.(J) 

/'  ^ 1+2ae^^' 


4N^(J) 


( J ) 


Equation  (87)^  for  all  points  2^  J is  solved  using  a Newton-Raphson 
iteration  method.  To  start  the  iteration  procedure^  a first  guess  for  ^(J) 
is  required  at  each  spatial  point.  This  guess  is  obtained  by  assuming 


charge  neutrality  throughout  the  whole  device.  Although  this  assumption  is 


not  valid  outside  the  bulk  regions,  the  solution  of  ♦(J)  in  the  transition 


regions  will  be  corrected  after  a few  iterations. 


Therefore  at  each  spatial  point  the  following  non-linear  equation  is 


solved  using  Newton's  Method. 


N.(J) 


FC4(J)3 


-a<J)  - -*(J)  . "O""  _ 

= n.e  ’ -n.e'’  + - 

' ^ 1+2ae*^'’^ 


4N^(J) 


n 


A first  guess  for  Newton's  method  is  easily  obtained  as  follows: 
For  N-type  regions,  Np(J)  * N(J)  and  ~ 0 


Figure  10,  Spaclal  steps  for  thermal 
equil Ibrium  case 
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N (J) 

♦ (J)  = ln(-^ ) 

"i 


(89a) 


Mhile  for  P-type  regions^  ***A^**^  “ ~ ® 


♦'J'  = 


(89b} 


Starting  with  either  equation  (89a),  or  (89b)  the  correction  in  ♦(J),  for 


each  iteration  is: 


A ♦ (J)  = -FC*„(J)]/F  C*„(J): 
K+1  ^ 


The  iteration  is  repeated  until  At|^(J)/4|^^^(J)  is  very  snail.  Then  having 
the  space  charge  neutrality  solution,  the  Newton-Raphson  method  is  applied 
to  solve  equation  (87)  as  follows: 


TftJTT 


AjA4(J1)  + BjA4(J)  + CjA4(J2)  = Dj  (91) 

This  procedure  is  repeated  for  (M-1)  spatial  points.  Using  the  boundary 
condition  of  = 0 at  the  ohmic  contact  or: 


A^d)  = A4(2) 


A4(N)  = A4(N-1) 


(92a) 


(92b) 


The  following  matrix  equation  results: 


^ ^ L—J  [v,J 

Equation  (93)  is  solved  for  each  iteration  by  using  the  algorithm  of  Appen- 
dix B.  The  iteration  is  continued  until  A^(J)/#(J)  reaches  a very  small 
number  for  2 £ J £ M-1.  Knowing  the  exact  spatial  distribution  of  ♦(x),  the 
electron  and  hole  distributions  are  easily  obtained  from  equations  (79)  and 
(80). 

To  calculate  the  electric  field,  one  has  to  keep  in  mind  that  the  zero 
electron  and  hole  currents  are  the  main  numerical  constraint  for  thermal 
equilibrium.  At  this  point  we  have  to  make  a decision  on  whether  to  use 
different  space-charge  points  and  field  points,  or  to  consider  all  points  as 
space  charge  and  field  points.  Also  a decision  has  to  be  made  on  whether  to 
calculate  the  electric  field  directly  from  Poisson's  equation,  or  to  calcu- 
late it  from  = 0 and  Jp  = 0.  In  the  following  discussion,  these  four 
possibilities  are  investigated  and  compared  with  each  other. 

(a)  Using  the  same  points  as  space  charge  points  and  f iel d points; 


1 
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A4(J1)  - ♦(J)  - ♦(JD 


A«(J)  = *(J2)  - 4(J) 


3N(J)  _ N(J2)  - N(J1)  ^ "i  r ♦(J)  + A* 
3x  2ax  2a X L 


(J)  _ g^CJ)  - A+(J1) 


- w(j)  r,A*(j)  _ -A* (Jin 

- TET  L'  ' J 


(94) 


3P(J)  I’g-^CJ)  - A4(J)  _ g-4(J)  + A4(Jl)‘j 


3x  2ax 


J 


_ P(J)  j■g-A♦(J)  _ gA4(Jl)J 


(95) 


JN(J)  = qvi^CE(J)«N(J)  + 


(96) 


VT«P(J)  ,.AA(J1)  .-A4(J),^ 

JP(J)  = qUpCE(J)*P(J)  + — 2a^ — -e  )J 


(97) 


(i)  Using  Poisson's  Equation; 


E(J) 


= IIl  = yT(AA(Jl)  •»  A»(J)) 


3x 


TKT 


(98) 


If  A4(J)  and  A4(J1)  « the  exponential  terns  of  equations  (96)  and  (97) 
can  be  approxinated  by  a linear  tern  and  only  then  JN  and  JP  will  have  a 
zero  value. 

The  total  potential  variation  across  a typical  transition  region,  as 
shown  in  Fig.  8,  is  about  0.8  volts.  If  this  region  is  divided  into  a 100 
spatial  points,  the  average  value  of  A4(J)  will  be  given  by: 


•r 
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A4(J)  x = 0.3  which  is  not  « 1 


Therefore^  this  method  requires  a very  large  number  of  spatial  points^  which 
makes  it  difficult  to  use. 


(ii)  Using  the  current  equations: 


In  order  to  calculate  the  electric  field  from  equations  (96)  and 
(97)  so  that  JN  - JP  s 0/  A4(J)  has  to  be  equal  to  A4(J1)^  which  is  not  the 
case,  expecially  if  a non-uniform  mesh  is  used. 


(b)  Different  space-charge  points  and  f iel d points: 

Following  the  same  argument  as  in  (a)  and  using  Fig.  12,  we  get: 


3N(1)  . ^ re*^**^  * 
ax  ■ Ax  ^ ® 


_ N(J)  r-AA(J) 


aP(I)  _ i r -♦(J2)  _ -♦(J2)  + A4(J)t 
ax  ■ Ax  '■  ■' 


P(J2)  ..  _ ^♦(j). 
Ax  ■* 


(100) 


for  E < 0: 


JN(I)  = 


qu  CE(I)«N(J)  ♦ lElfliJl  (e^*^**^-1)] 
n AX 


JP(I)  = qiipCE(I)*P(J2)  ♦ (e^*^’*^-1)D 


(i)  Using  Poisson's  Equation: 


(101) 


(102) 
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Ed)  = 


(103) 


As  in  case  (a-i),  the  use  of  Poisson's  equation  requires  a very  large 
nunber  of  spatial  points  to  insure  that  JN  > JP  = 0. 


(ii)  Using  the  current  equations: 

From  equations  (101)  and  (102)/  one  can  easily  obtain  a value  for 
the  electric  field  to  satisfy  the  zero  current  requirement 


Ed)  = (1  - 

AX 


(104) 


Therefore  the  electric  field  is  calculated  using  equation  (104). 

It  should  be  noted  that  for  snail  values  of  44(J)/  equation  (104)  can 
be  approximated  by  equation  (103).  Therefore/  in  order  to  satisfy  both  the 
zero  current  condition  and  Poisson's  equation/  4^(J)  should  have  a small 
value  inside  the  P*-N  transition  region.  As  a result/  a uniform  space-mesh 
is  safe  to  use  in  the  case  of  a low  doped  W-region/  while  for  high  doping 
where  the  potential  variation  in  the  P^  -N  transition  region  is  large,  a 
nonuniform  space-mesh  is  necessary. 

Although  all  the  formulations  and  numerical  methods  apply  for  any  two 
terminal  one-dimensional  semiconductor  device/  from  now  on  we  will  be  con- 
cerned with  discussing  and  presenting  the  results  for  the  P^-N-N*  diode 
structure. 


Device  under  consideration; 

The  typical  construction  of  a semiconductor  P*-N-N*  diode  is  illustrat- 
ed in  Fig.  13.  The  junction  may  be  circular  with  a diameter  of  3 mils  to  5 
mils.  The  lightly  doped  N“  bulk  semiconductor  may  be  doped  from  10^^/cm^  to 
10^^/cm^  and  has  a thickness  (w)  from  10  microns  to  300  microns.  The 
P*-region  is  heavily  doped  from  10^®/cm^  to-  10^^/cm^.  The  N*-region  is 


Large  Area  Ohmic  Contact 


Figure  13 


Cross  section  of  a typical  P N P junction 
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doped  at  about  10  tcnr,  and  its  thickness  (b)  may  vary  between  10  and  25 
mils. 

The  doping  profile  of  both  the  p'*'  and  the  N*  regions  is  approximated  by 
the  complementary  error  function  distribution  as  in  equation  (105). 

C(x,t)  = C.  erfc  ( - -)  (105) 

* 2V0t 

where: 

Cg  is  the  constant  surface  concentration 

D is  the  diffusion  coefficient 

t is  the  diffusion  time  required  for  C(x^t)  to 

reach  the  doping  concentration  of  the  N-type 

at  the  physical  junction  location  x. 

In  this  chapter^  the  thermal  equilibrium  solution  is  presented  for  two 

14  3 

device  structures.  The  first  device  has  a low  doped  N region  of  10  /cm 
and  a total  length  of  10  microns.  The  metallurgical  junction  of  the  P^-N 
transition  region  is  located  at  a distance  of  1.6  microns.  Figure  14  shows 
the  uniform  space  mesh  used  with  198  spatial  steps,  198  field  points,  and 
199  space  charge  points;  only  the  space  charge  points  are  indicated  in  Fig. 
14.  Figures  15  and  16  show  the  electric  field  spatial  distribution  and  the 
carrier  density  distributions. 

The  second  device  has  a total  length  of  20  microns.  The  doping  density 
of  the  N-region  is  chosen  to  be  equal  to  2.64  x 10^^/cm^  in  order  to  match 
the  fabricated  device  discussed  in  Chapter  VI.  A non-uniform  space-mesh  is 
used  with  239  spatial  steps,  239  field  points,  and  240  space  charge  points. 
As  discussed  before,  the  point  density  is  highest  inside  the  P*-N  transition 
region,  to  insure  a small  potential  variation  "A4^(J)".  Figures  17  and  18 
show  the  electric  field  distribution  and  the  carrier  density  distributions 
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Figure  16.  Electron  and  hole  density  distribution  of  thermal 
equilibrium  for  diode  S\ 
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CHAPTER  IV 
REVERSE  PULSED  CASE 

In  this  chapter,  the  numerical  technique  and  the  formulations  of 

Chapter  II  are  used  to  solve  for  the  transient  behavior  of  the  two  P*-N-N* 

structures  previously  considered  in  Chapter  III. 

14  3 

The  device  of  10  microns  length  and  10  /cm  doping  density  is  con- 
sidered first.  The  thermal  equilibrium  conditions  of  Figs.  15  and  16  in 

Chapter  III  are  used  here  as  the  initial  conditions.  The  external  excita- 

2 

tion  is  in  the  form  of  a negative  current  step,  of  -10  amp/cm  amplitude, 
and  0.1  nanosecond  rise  time.  The  numerical  solution  simulates  the  device 
transient  behavior,  between  the  thermal  equilibrium  and  the  steady  state. 
The  time  step  "A  t"  is  chosen,  when  implementing  the  computer  program,  in 
such  a way  that  a converged  solution  is  obtained  after  a maximum  of  five 
iterations  for  each  instant  of  time.  If  not,  a new  At  is  chosen. 

Fig.  19  shows  the  applied  current  density  vs.  time,  along  with  the 
ouput  voltage  vs.  time.  The  time  axis  is  divided  into  three  regions.  A,  B, 
and  C.  Region  A describes  the  behavior  between  thermal  equilibrium  and 
punch  through  of  the  N”-region.  Region  B is  the  period  between  punch 
through  and  the  beginning  of  impact  ionization.  Region  C describes  the 
avalanching  period  uniti  a steady  state  solution  is  reached. 

The  spatial  distributions  of  the  electrons  and  holes,  displacement 
current  density,  carrier  current  densities,  and  electric  field  for  the  in- 
stants of  time  indicated  in  Fig.  19  are  shown  in  Figs.  20  through  38. 

Til#  transient  solution  of  region  A is  shown  in  Figs.  20  to  23.  Immedi- 
ate* r after  the  current  step  application,  the  electric  field  will  change  in 
a aariod  of  time  in  order  to  support  the  total  current.  Therefore, 
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the  displacement  current  component  will  be  equal  to  the  total  current 
throughout  the  whole  device  according  to  the  following  equation: 


J(T=0^  = c|| 


Shortly  after  the  step  application,  the  build-up  of  ohmic  voltage  across  the 
heavily  doped  P*-region  causes  the  total  current  to  be  carried  by  the  major- 
ity carriers  (holes)  instead  of  being  carried  by  the  displacement  current. 
Similarly,  in  the  N~-N*  region,  the  total  current  will  be  carried  by  majori- 
ty carriers  (electrons).  In  the  P^-N  depletion  region,  the  carrier  densi- 
ties are  very  small,  too  small  to  support  any  current.  Therefore,  the  total 


current  will  still  be  carried  by  the  displacement  current  component  as  shown 
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in  Fig.  21.  At  this  stage,  the 


3x 


term  of  the  electron  continuity 


equation  near  the  depletion  region  is  negative  in  the  N - region.  Also,  the 


1 
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3x 


term  of  the  hole  continuity  equation  is  negative  between  the  P - 
region  and  the  depletion  region  as  shown  in  Fig.  20  for  times  t^  to  t^. 
Mathematically  this  causes  recombination  of  majority  carriers  near  the  de- 
pletion region.  Physically  it  represents  extraction  of  majority  carriers, 
or  in  other  words,  depletion.  This  carrier  depletion  is  shown  in  Fig.  22 
for  region  A.  As  the  N”-region  continues  to  deplete,  the  displacement 
current  component  replaces  the  electron  current  component.  Therefore,  the 
electric  field  increases  negatively  in  the  depletion  region  until  punch- 
through  at  time  t6.  These  electric  field  distributions  are  shown  in  Fig.  23 
for  region  A. 

In  Region  B,  the  current  is  mainly  carried  by  the  displacement  current 
in  the  N -depleted  region  as  shown  in  Fig.  25.  The  carrier  depletion  con- 
tinues at  a much  slower  rate  as  shown  in  Fig.  25.  The  electric  field  con- 
tinues to  increase  negatively  until  it  approaches  the  avalanching  value  of 
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about  -2x10^  volt/ca  at  tiae  t^j  as  shown  in  Fig.  27. 

Region  C describes  the  avalanching  part  of  the  transient  behavior  as 
explained  graphically  in  Figs.  28  to  35.  In  the  first  part  of  this  region, 
the  electron  and  hole  impact  ionization  coefficients  increase  very  rapidly 
as  the  electric  field  continues  to  increase. 

At  this  point,  it  should  be  noted  that  since  no  acceleration  terms  are 
present  in  the  equations,  the  response  of  the  carriers  to  the  electric  field' 
in  a computer  simulation  is  instantaneous.  Whereas,  it  is  known  that  some 
time  delay  does  occur  due  to  the  fact  that  a particle  needs  time  to  ac- 
celerate and  fully  respond  to  the  applied  field. 

Therefore,  the  rapid  increase  of  the  impact  ionization  coefficients 
causes  the  carrier  densities  to  reach  a high  level.  These  large  carrier 
densities  combined  with  the  high  electric  field  result  in  the  particle 
current  density  components  assuming  values  which  are  higher  than  the  exter- 
nally applied  current  density.  This  in  turn  causes  the  displacement  current 
to  switch  to  a positive  value.  This  situation  is  illustrated  by  the  curves 

of  the  time  instant  t^^  in  Figs.  28  to  31. 

As  soon  as  the  magnitude  of  the  electric  field  drops,  the  impact  ioni- 
zation coefficients  and  the  particle  currents  decrease,  and  again  a large 
portion  of  the  total  current  is  carried  by  the  displacement  current  com- 
ponent as  illustrated  by  the  curves  of  the  time  instant  t2Q  Figs.  32  to 
35.  This  trade-off  between  the  particle  currents  and  the  displacement 
current  causes  the  output  voltage  to  fluctuate  around  an  average  value  as 
shown  in  Fig.  19  until  near  steady  state  where  the  particle  current  densi- 
ties replace  the  displacement  current  component  as  shown  in  Figs.  32  and  33. 
At  steady  state,  the  displacement  current  component  vanishes  throughout  the 
whole  device.  Also,  the  generation  rate  due  to  the  avalanching  mechanism  is 
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Figure  25.  Displacement  current  density  distributions  for  time 
Instants  t_  to  t,,  as  shown  In  Figure  19 
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Electron  and  hole  density  distributions  at  time  instants 
t-  to  t,,  as  shown  in  Figure  19 
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balanced  by  the  rate  of  carriers  extraction,  and  the  ShockI ey-Read-Hall 
recombination  rate.  In  other  words,  the  two  continuity  equations  reduce  to 
the  following  two  equations. 


{ I 


i 

® i TIT 


(106) 


G = -R  - ^ 


q ax 


(107) 


where; 


G = aUnl  + BiJpl 


pn  - n. 
Tg(p+n+2n^.) 


From  equations  (106)  and  (107),  it  is  clear  that  at  steady  state  the 
following  condition  has  to  be  satisfied 
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The  steady  state  solution  is  shown  in  Figs.  36  to  38.  The  solution  in 
Fig.  19  is  obtained  by  dividing  the  time  axis  into  a thousand  non-uniformi y 
spaced  time  steps.  The  time  steps  have  relatively  small  values  in  region  A, 
where  the  time  variations  of  al I the  physical  quantities  are  changing  very 
rapidly.  In  region  B,  the  time  steps  are  made  larger  since  the  electric 
field  is  mainly  the  variable  quantity.  In  region  C,  the  time  step  is  re- 
duced again  in  the  avalanching  period.  As  the  steady  state  approaches,  the 
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Figure  32,  Electron  and  hole  current  density  distributions  for  time 


Instants  t-Q  to  t22  as  shown  in  Figure  19 
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tiae  step  could  be  adjusted  to  a very  large  value  without  any  loss  of  accu- 
racy. This  is  one  of  the  very  iaportant  advantages  of  using  the  iterative 
iaplicit  technique. 

The  tiae  levels  of  the  different  points  indicated  in  Fig.  19,  along 
with  their  tiae  step  nuaber  and  the  output  voltage  are  suamarized  in  Table 
2. 

Siailarly^  the  transient  behavior  of  the  second  device  of  20  ai crons 

length  and  2.64  x 10^^/ca^  doping  density  is  presented  in  Figs.  39  to  51. 

The  theraal  equilibriua  conditions  of  Figs.  17  and  18  in  Chapter  III  are 

used  here  to  initiate  the  transient  solution.  The  external  excitation  is  in 

2 

the  fora  of  a negative  current  step  of  -100  aap/ca  aapi itude,  and  0.1  nano 
second  rise  tiae. 

Figue  39  shows  the  external  applied  current  density  vs.  tiae,  along 
with  the  output  voltage  vs.  tiae.  The  spatial  distributions  of  the  electron 
and  hole  current  densities,  the  displacement  current  density,  the  carrier 
densities,  and  the  electric  field  for  the  instants  of  time  Indicated  in  Fig.  | 

39  are  shown  in  Figs.  40  to  51.  There  are  mainly  three  differences  between 
this  device  and  the  first  device  previously  considered: 

(1)  longer  device 

(2)  higher  doping  in  the  N -region 

(3)  larger  amplitude  of  the  externally  applied  current  step. 

The  second  difference  causes  a slower  rate  of  the  N -region  depletion, 
which  is  combined  with  the  first  difference  result  in  a partial  depletion  of 
the  N”-region  at  steady  state.  The  second  and  third  differences  are  respon- 
sible for  the  higher  electric  field,  and  for  more  voltage  fluctuation  near 
steady  state.  I 
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Figure  3^».  Electron  and  hole  density  distributions  for  time 
Instants  t2Q  to  t22  as  shown  in  Figure  19 
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Figure  36,  Electron  and  hole  current  density  distributions  for  time 


Instant  t^,  as  siown  In  Figure  19 
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Figure  37.  Electron  and  hole  density  distributions  for  time 
Instant  t.,  as  shown  In  Figure  19 


- 94  - 


Figures  40  to  43  demonstrate  a very  interesting  phenomena  in  the  sense 
that  after  time  t^  (=1.8  nanoseconds)  both  carrier  depletion  and  avalanche 
multiplication  are  simultaneously  present.  This  is  due  to  the  triangular 
shape  of  the  electric  field.  Therefore^  the  high  electric  field  near  the 
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tion^  cause  carrier  depletion.  This  is  very  clearly  illustrated  in  the 
curves  of  the  time  instant  tj.  Figs.  44  to  47  show  the  spatial  distribu- 
tions of  the  different  quantities  at  one  of  the  voltage  fluctuating  time  in- 
stants t^  (=  4.8  nanoseconds)^  while  Figs.  48  to  51  show  the  steady  state 
solution  at  time  instant  t^  (=8  nanoseconds).  The  time  levels  of  the  dif- 
ferent points  indicated  in  Fig.  39/  along  with  their  time  step  number  and 
the  output  voltage/  are  summarized  in  Table  3. 

To  summarize  the  reverse  pulsed  conditions/  the  device  with  low  doping 
density  and  short  N -region  depletes  at  a relatively  fast  rate  until  punch 
through  is  reached.  The  electric  field  increases  uniformly  in  the 
N -region/  until  the  onset  of  avalanche  by  impact  ionization.  On  the  other 
hand/  the  device  with  high  doping  density  and  long  N~-region  does  not  fuMy 
deplete.  In  this  case  the  electric  field  increases  while  keeping  its  tri- 
angular shape  in  the  depletion  region  until  it  reaches  avalanche.  In  both 
caseS/  the  steady  state  condition  is  satisfied  when  the  rate  of  generated 
carriers  via  avalanche  multiplication  is  balanced  by  the  recombination  and 
depletion  rates  as  shown  in  equations  (106)  and  (107). 


Time  In  nanoseconds 

Figure  39.  Current  density  and  output  voltage  for  the  reverse  biased  case  (second  device) 
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Figure  40,  Electron  and  hole  current  density  distributions  for 
time  instants  tj  to  t^  as  shown  in  Figure  39 
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Figure  43.  Electric  field  distribution  for  time  instants  t 
as  shown  In  Figure  39 
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Figure  48. 


Electron  and  hole  current  density  d 
time  instant  t^  as  shown  in  Figure  3l 
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Figure  50 


Electric  field  distribution  for  time  Instant  t 
In  Figure  39 


as  shown' 
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Table  3«  Output  voltage  as  a function  of  time  for  Figure  39 


Time  Level 


Time  in 

Nano  Seconds 

Number  of 

Time  Steps 

Output  Voltage 

In  Vol ts 

0.6 

60 

- 7.913 

1.2 

120 

- 22.895 

1.8 

180 

- 45.16 

3.0 

2k0 

-117.86 

3.6 

300 

-166.68 

l».8 

420 

-142.99 

8.0 

560 

-150.2 
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Figure  51.  Pure  implicit  technique  flow  chart 
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CHAPTER  V 

FORWARD  PULSED  CASE 

The  Main  purpose  of  this  chapter  is  to  explain  the  unexpectedly  high 
voltage  which  appears  across  the  P*-n"-N*  structure  when  pulsed  with  extreme 
current  densities  in  the  forward  pulsed  direction.  A qualitative  discussion 
is  presented  first,  followed  by  the  computer  simulation  results  using  the 
conventional  Shocki ey-Read-Hal I recombination  term  in  the  two  continuity 
equations.  Discussion  of  the  results  leads  to  the  prediction  of  an  effec- 
tive transient  lifetime  which  varies  with  time  and  reduces  to  its  equilibri- 
um value  near  the  steady  state  condition. 

In  the  P*-n”-N^  structure,  there  are  four  regions  which  may  give  rise 
to  the  unexpected  high  voltage  across  the  device,  the  heavily  doped  P^  and 
N*  regions,  the  P*-n"  transition  region,  and  the  n"  region.  The  heavily 
doped  regions  operate  under  low  level  injection  conditions;  therefore,  their 
ohmic  voltage  drop  is  almost  nonexistent.  The  voltage  across  the  P*-N 
transition  layer  is  very  small  due  to  the  forward  biased  condition.  There- 
fore, the  low-doped  N^region  must  be  considered  the  cause  of  the  large  vol- 
tage. This  can  be  explained  by  considering  the  behavior  of  an  N-type  sem- 
iconductor at  an  extreme  current  density.  At  relatively  low  current  densi- 
ties bulk  silicon  will  obey  Ohm's  Law  and  the  current  density  is  given  by 
equation  (109)  and  (110). 


(109) 


1 


J = qNv^ 


(110) 


where, 

N = electron  concentration,  #/cin^ 

-19 

q = 1.602  x 10  , coulomb 

2 

11^  = mobility,  cm  /v-sec 

t = electric  field,  v/cm 

= average  drift  velocity  of  the  carriers,  cm/sec 

This  relationship  is  valid  at  low  values  of  electric  field.  At  larger  elec- 
tric fields  (or  current  density)  the  mobility  will  become  field  dependent 
E5J  as  shown  in  Fig.  52. 

At  electric  fields  below  the  critical  value  for  avalanche  breakdown, 
the  electron  concentration  is  constant.  The  resistance  of  the  bulk  region, 
therefore,  increases  as  the  electric  field  increases.  Since  Fig.  52  is  for 
very  pure  material,  for  doped  material  one  would  expect  the  slope  to  be 
smaller  at  low  values  of  the  electric  field.  However,  in  both  cases  at  an 

4 

electric  field  of  about  2 x 10  v/cm  the  electron  velocity  approaches  a 
scattering-limited  value  of  about  10^  cm/sec.  Much  higher  electric  fields 
are  required  for  holes  in  silicon,  2 x 10^  v/cm  before  the  scattering  limit- 
ed velocity  is  obtained.  For  N-type  silicon  an  applied  field  in  excess  of 
10^  v/cm,  the  velocity  being  scattering-limited,  implies  a large  resistivi- 
ty. Or,  looking  at  it  from  the  current  standpoint,  by  using  equation  (110), 
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3 = qNVj  = 10^^  X 1.602  x 10“^^  x 10^  = 160  Amp/cm^ 

14  3 

in  the  high  resistance  region^  for  an  N-region  doped  at  10  /cm  , where 
N - below  avalanche.  This  implies  that  the  larger  the  N-region  is 
doped,  the  larger  the  current  density  needed  to  place  it  in  the  high  resis- 
tance region.  If  the  bulk  N-region  is  7 microns,  then, 

V = 2 X 10^  v/cm  X 7 X 10  ^ = 14  volts. 

This  voltage  appears  across  the  bulk  region  alone,  independent  of  the  junc- 
tion being  forward  or  reverse  biased. 

At  current  densities  larger  than  160  Amp/cm^,  v^  cannot  increase; 
therefore,  the  electric  field  increases  to  the  avalanching  value  to  generate 
enough  carriers  to  support  the  larger  current  density.  This  occurs  at  an 
electric  field  of  about  1.5  x 10^  V/cm  for  electrons  in  silicon.  In  this 
case,  the  voltage  across  the  bulk  region  is  at  least, 

V i 1.5  X 10^  X 7 X 10"^  = 105  volts. 


In  the  diode  structure,  these  bulk  effects  are  coupled  with  the  junction  ef- 
fects. As  discussed  in  Chapter  II,  the  time  variation  of  the  minority  car- 
ries follows  the  time  variation  of  the  electric  field  due  to  the  relaxation 
time  constant  being  shorter  than  the  transit  time. 

To  study  the  behavior  of  a P*-N-N^  diode  following  the  application  of  a 
low  level  current  pulse,  we  assume  that  the  diode  is  initially  at  thermal 
equilibrium.  To  insure  that  the  device  will  reach  a final  steady  state  be- 
fore the  pulse  end,  we  assume  that  the  pulse  width  is  greater  than  five 
times  the  transit  time  of  the  n"  region.  Immediately  after  the  pulse  appli- 
cation, the  electric  field  will  change  in  a short  period  of  time  in  ortter  to 
support  the  total  current.  Hence,  the  displacement  current  component  will 


be  equal  to  the  total  current  throughout  the  device  according  to  the  follow- 
ing equation: 

J(t  = O'^)  = e 41  (111) 


As  time  proceeds,  one  can  investigate  four  separate  regions.  After  a 
time  of  the  order  of  the  dielectric  relaxation  time  constant  in  the  two 
heavily  doped  regions,  the  build-up  of  ohmic  voltages  causes  the  total 
current  to  be  carried  by  the  majority  carriers  (holes  in  the  P*  region  and 
electrons  in  the  N*  region). 

Similarly,  in  the  low  doped  region,  the  majority  carrier  current  will 
replace  the  displacement  current  after  a period  of  time  of  the  order  of  the 
dielectric  relaxation  time  constant  in  the  N'-region.  In  the  P^-N  transi- 
tion layer,  the  minority  carriers  start  to  diffuse  and  charge  the  transition 
layer  capacitance.  This  phenomena  needs  a period  of  time  in  the  order  of 
the  transit  time,  during  which  the  hole  and  electron  currents  replace  the 
displacement  current  in  the  transition  layer. 

If  a large  current  pulse  is  applied  to  the  P*-N  -N^  diode,  initially 
the  voltage  that  would  appear  across  the  diode  will  be  the  sum  of  the  vol- 
tage across  the  transition  region  and  the  ohmic  voltage  drop  across  the  N 
region,  as  discussed  in  the  low  current  pulse  case,  except  that  the  ohmic 
drop  in  this  case  is  much  larger  due  to  lower  mobility  of  the  carriers  at 
the  higher  electric  field.  As  the  minority  carriers  start  to  diffuse,  the 
N**  region  will  operate  under  high  injection  conditions,  conductivity  modula- 
tion will  start,  and  the  voltage  across  the  N“-region  will  begin  to  decay 
with  a time  constant  approximately  equal  to  the  transit  time. 

If  an  extremely  intense  current  pulse  is  applied,  the  large  current 
density  will  initially  be  carried  by  the  displacement  current  component. 
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The  majority  carriers  in  the  N“  region  will  not  be  sufficient  to  replace  the 
displacement  current;  therefore,  the  electric  field  increases  very  rapidly 
to  create  avalanche,  which  in  turn  generates  enough  majority  and  minority 
carriers  to  replace  the  displacement  current.  In  this  case  the  N region  is 
under  high  level  injection,  and  yet  the  voltage  appearing  across  the  device 
is  very  high.  As  the  injection  and  diffusion  process  has  time  to  take 

place,  the  electric  field  necessary  to  support  the  current  decreases. 

In  the  remainder  of  this  chapter  the  numerical  method  of  Chapter  II  is 
used  to  solve  for  the  transient  behavior  of  the  p'*’-N-N^  structure,  with 

external  excitations  in  the  form  of  current  steps  with  different  amplitudes 
and  rise  times. 

The  device  under  consideration  is  the  first  device  presented  in  Chapter 

14  3 

III  with  a low  doping  region  of  10  /cm  , and  7 microns  width.  The  total 
device  length  is  10  microns.  The  N -region  "lifetime"  is  4 x 10  seconds. 

External  excitation  in  the  form  of  £ I ow  current  step 

A relatively  low  current  step  is  first  applied  to  the  device  initially 
under  thermal  equilibrium.  The  thermal  equilibrium  solution  is  described  by 
Figs.  15  and  16  in  Chapter  III.  The  current  step  has  a rise  time  of  1.25 

nanoseconds,  and  an  ampi itude  of  1 ampere/cm  . The  response  of  the  device 

is  illustrated  in  Figs.  53  to  57.  The  terminal  voltage  v(t)  as  a function 
of  time  is  shown  in  Fig.  53  along  with  the  external  excitation  J(t).  Figs. 
54  to  57  show  the  electron  and  hole  current  density  components,  the  dis- 
placement current  density  component,  the  carrier  densities,  and  the  electric 
field  as  functions  of  position  for  various  instants  of  time. 

Immediately  after  the  current  step  application,  the  current  will  be  to- 
tally carried  by  the  displacement  current  component  according  to  equation 
(111).  As  soon  as  any  electric  field  is  established  in  the  heavily  doped 


- 115  - 


e 

i 

I 

i' 

! 

r 

I 

I 

<• 

. i 

f 

t I 

> 

t 

. 5 

; i 

i 

' i 

^ I 

\ 

' i 

( 

: > 


regions,  the  current  will  be  carried  by  the  majority  current  component. 
This  situation  is  illustrated  in  Figs.  54  and  55  for  the  time  instant  t^ 
where  the  total  current  is  equal  to  0.2  amp/cm^.  At  this  time  level,  the 
current  is  partially  carried  by  the  majority  carrier  component,  and  partial- 
ly by  the  displacement  current  component  in  the  low  doped  region,  while  it 
is  totally  carried  by  the  displacement  current  component  in  the  depletion 
region.  A positive  displacement  current  means  an  increase  in  the  electric 
field,  or  less  negative  value  as  shown  in  Fig.  57.  After  a short  time,  the 
current  in  the  N~-region  will  be  totally  carried  by  the  majority  carrier 
component.  This  occurs  within  a period  approximately  equal  to  the  dielec- 
tric relaxation  constant  of  the  N“-region.  This  is  illustrated  by  the 
curves  of  the  time  instant  where  the  displacement  current  is  still  equal 
to  the  total  current  in  the  depletion  region.  At  time  level  t^,  the  deple- 
tion region  becomes  narrower  due  to  the  increase  in  the  electric  field,  and 
to  some  carrier  diffusion  across  the  junction.  This  causes  the  electron 
current  component  to  further  replace  the  displacement  current  component  on 
the  right  side  of  the  depletion  region.  At  this  point,  it  should  be  noted 
that  one  of  the  main  differences  between  a forward  biased  junction  and  a re- 
verse biased  one,  is  the  sign  of  the  spatial  derivatives  of  the  electron  and 
hole  currents  near  the  depletion  region.  In  Chapter  IV  it  was  shown  that 
these  two  terms  cause  carrier  recombination  or  in  other  words  carrier  deple- 
tion near  the  reverse  biased  junction.  From  Fig.  54,  it  is  clear  that  the 

JL  liill  and  -—  terms  of  the  particle  continuity  equations  are  posi- 
q |3x  I q |3x  I 

tive  near  the  depletion  region.  This  causes  electron  diffusion  from  the 
N”-region  to  the  P^-region,  and  hole  diffusion  from  the  P^-region  to  the 
N”-region.  As  time  proceeds,  these  carrier  diffusions  increase  and  cause 
the  electric  field  to  decrease  in  the  N"-region  and  also  cause  the  displace- 
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Figure  53.  J(t)  and  V(t)  for  a rela 
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merit  current  to  decrease  in  the  depletion  region.  This  process  continues 
until  steady  state,  where  the  displacement  current  goes  to  zero  throughout 
the  entire  device,  and  the  total  current  is  carried  by  both  electrons  and 
holes  as  illustrated  by  the  curves  of  the  time  instant  t^  (which  is  near 
steady  state).  It  should  be  noticed  that  after  a time  of  the  order  of  t2, 
the  displacement  current  in  the  N-region  is  almost  non-existent  compared  to 
that  in  the  depletion  region.  Therefore,  any  time  variation  of  electric 
field  inside  the  N -region  does  not  affect  the  output  voltage  for  this  level 
of  injection. 

External  excitation  in  the  form  of  extremel y high  current  step 

4 2 

An  extremely  high  current  step  of  10  Amp/ cm  amplitude  and  two  pico 
seconds  rise  time  is  applied  to  the  P^-N  -N^  structure.  The  computer  simu- 
lation results  are  illustrated  in  Figs.  58  to  69  for  nine  instants  of  time. 
Shortly  after  the  step  application  (t^  and  t2)^  the  current  in  the  heavily 
doped  regions  is  carried  by  the  majority  carrier  components,  while  the  dis- 
placement current  equals  the  total  current  in  the  low  doped  and  the  depi e- 

1 I d jpI 

tion  regions.  As  shown  in  Fig.  59,  both  quantities -r — near  the 

Q I d X I 

near  the  N*-region  are  positive  and  have  very  large 
values  due  to  the  high  current  level.  On  the  other  hand,  the  quantities 

^ 1 I d jnI  ^ 

near  the  N -region  and  — l-r — near  the  P -region  are  very  smal  I . 

q 1 3x  I 

This  causes  a very  large  number  of  holes  to  diffuse  from  the  P -region  to 
the  N^-region,  and  a similar  amount  of  electrons  to  diffuse  from  the 
N*-region  to  the  N”-region.  In  other  words,  the  P*-N  and  the  N*-N  transi- 
tion regions  show  identical  behavior.  This  increase  of  the  diffused  car- 
riers causes  the  hole  current  to  replace  the  displacement  current  in  the 
left  side  of  the  N’-region,  while  the  electron  current  replaces  it  in  the 
right  side.  Also  it  is  clear  from  Fig.  62  that  the  electric  field  increases 
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P -region,  and 
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very  rapidly  with  time,  due  to  the  large  displacenent  current  value,  until 
it  reaches  an  avalanching  value.  At  first  (t^  and  t2)  the  avalanche  multi- 
plication effect  is  noticeable  in  the  hole  increment.  As  the  electric 
field,  the  electron  current,  and  the  hole  current  increase,  the  avalanche 
multiplication  gets  larger.  As  a result,  the  avalanched  carriers  start  to 
play  a role  in  current  transport  as  shown  for  time  instant  tj.  This  process 
continues  in  Figs.  62  and  64.  The  carrier  levels  increase  due  to  diffusion 
and  avalanche  multiplication,  the  particle  currents  increase,  and  as  a 
result,  the  displacement  current  decreases,  then  switches  to  a negative 
value  as  the  particle  currents  become  larger  than  the  external  applied 
current.  This  causes  the  electric  field  to  decrease  in  the  N -region,  and 
consequently,  the  output  voltage  decreases  as  shown  in  Fig.  S8. 

To  reach  a steady  state,  the  displacement  current  and  the  time  varia- 
tion of  both  holes  and  electrons  should  go  to  zero  values  or, 

JN  + JP  = Jt  (112) 


PN  - n? 
Tg(P+N+2n^) 
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_ 1 

ajN 
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ax 

“ q 

ax 

(113) 


equation  (112)  is  satisfied  for  time  instants  tg  and  t^  as  shown  in  Fig.  67, 
but  the  condition  of  the  following  equation  is  also  satisfied. 


1 lliil  = I ~ ^i 

q lax  I ~ q lax  I TQ(P+N+2n^) 


(114) 


This  is  again  due  to  the  high  current  values.  Equation  (114)  causes 
the  carrier  levels  to  keep  rising  throughout  the  entire  device  until  the 
minority  carriers  in  the  heavily  doped  regions  become  comparable  to  the  dop- 
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Figure  62,  Electric  field  distribution  for  time  Instants  t = 0 to  t 
as  shown  In  Figure  58 
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Figure  63 


Electron  and  hole  current  density  distributions  for  time 
Instants  t^  to  t^  as  shown  In  Figure  58 
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ing  density,  which  causes  the  current  at  the  ohmic  contacts  to  be  carried  by 

both  the  minority  and  the  majority  carriers,  in  such  a way  as  to  reduce  the 

} 

particle  current  derivatives  and  to  satisfy  equation  (113).  This  situation 
continues  until  t^  where  the  output  voltage  drops  to  a value  of  2.1  volts 
without  reaching  a steady  state. 

From  the  results  presented  previously,  it  is  clear  that  although  the 
present  formulations  are  adequate  in  explaining  the  transient  response  for 
relatively  low  current  densities,  it  fails  to  give  a physical  explanation 
for  extreme  current  levels,  where  the  device  is  too  far  from  the  equilibrium 
condition.  This  is  mainly  because  in  the  forward  bias  case,  the  present 
formulation  does  not  provide  a means  of  extracting  any  plasma  formation  in 
the  N”-region.  Having  that  in  mind,  the  transient  response  of  the  previous 
section  is  reanalysed  during  three  stages,  the  excess  diffusion  stage  near 
t=Q,  the  impact  ionization  stage,  and  near  steady  state  where  the  displace- 
ment current  is  very  small. 


The  excess  diffusion  stage  and  the  transient  I ifetime  prediction 

The  transient  behavior  near  t = 0,  for  example  time  instants  t^  and  t^, 
is  shown  in  Figs.  59  to  62.  The  time  rates  controlling  the  electron  and 
hole  changes  are  described  by  the  two  continuity  equations: 
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where 


and  G assumes  a small  value  at  this  stage. 


The  extremely  high  value  of  compared  to  in  the  P*-n"  tran- 
sition region  causes  the  number  of  holes  to  increase  at  a much  higher  rate 
compared  to  the  increase  in  electrons.  A similar  situation  develops  in  the 
N*-N~  transition  region  causing  a large  increase  of  electrons  compared  to 
holes.  Whereas,  it  is  known  that  in  a neutral  or  in  a quasi-neutral  region, 
there  is  the  restraint  on  the  system  that  the  net  change  in  negative  charge 
is  balanced  by  a change  in  the  positive  charge  C7J.  It  is  also  known  that 
in  case  of  low  level  injection,  due  to  the  complicated  nature  of  the  two 
continuity  equations,  it  is  customary  to  solve  a simplified  continuity  equa- 
tion which  controls  the  time  variation  of  the  minority  carriers.  Then  the 
assumption  of  a N = A P is  usually  used. 

The  above  discussion  suggests  that  the  two  recombination  terms  in  equa- 
tions (115)  and  (116)  may  not  be  equal,  especially  when  the  device  is  too 
far  away  from  equilibrium,  as  in  the  present  case.  Blakemore  C8]  pointed 
out  that  these  two  recombination  terms  can  be  different  in  the  two  continui- 
ty equations,  and  that  there  is  a hope  of  finding  equal  terms  only  in  care- 
fully purified  semiconductors  of  small  intrinsic  band  gap  width  at  high  tem- 
peratures, when  band-to-band  processes  dominate  the  recombination  behavior. 
Therefore,  the  two  continuity  equations  are  rewritten  with  different  recom- 
bination terms  as  follows: 
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T^(P+N+2n,) 


The  transient  lifetimes  t and  -t  are  adjusted  in  such  a way  to  limit 


the  increase  of  majority  cas  to  that  of  minority  carriers  in  the  two  heavily 
doped  regions,  which  normally  operate  under  low  level  injection.  The  heavi- 
ly doped  regions  are  defined  according  to  the  difference  between  the  two 


near  the  ohmic  contacts  as  follows 


quantities 


P -region  where 


(119) 


N -region  where 


Throughout  the  entire  device  the  minority  carrier  continuity  equation 


assumes  the  normal  lifetime  "t 


while  the  transient  lifetime  in  the 


ity  carrier  continuity  equation  is  adjusted  to  satisfy  equation  (121)  in 


regions 


Therefore,  in  the  heavily  doped  regions,  equations  (117)  and  (118)  are 


equated  to  give  the  following  transient  lifetimes 
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and,  in  the  N -region. 
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where 


PN  - n^ 
P+N+2n7 


Tg  = equilibrium  lifetime. 

The  lifetime  variation  inside  the  device  is  shown  in  Fig.  70.  At  the 

start  of  the  transient  solution,  the  transient  lifetime  is  very  short,  and 

Id  JPI 
Tr~\ 

o X I 

approaches  transient  lifetime  reduces  to  the  equilibrium 

lifetime  "tq"  throughout  the  entire  device. 

The  impact  ionization  stage 

As  soon  as  the  electric  field  reaches  a value  of  about  2 x 10^  volt/cm, 
the  electron  and  hole  levels  start  to  increase  in  the  N -region,  which 
causes  an  increase  in  electron  and  hole  currents.  The  particle  currents  in- 
crease, combined  with  the  electric  field  increase  result  in  a higher  rate  of 
carrier  generation.  Since  in  the  present  formulation  there  is  no  means  to 
extract  the  excess  carriers,  the  electron  and  hole  currents  continue  to  in- 
crease until  they  exceed  the  external  applied  current.  This  causes  the  dis- 
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placement  current  to  switch  to  a negative  value;  and  while  the  electric 
field  is  d'creasing^  the  carrier  level  keeps  increasing  until  the  electric 
field  drops  below  the  avalanching  value  of  2 x 10^  volt/cm.  Therefore,  a 
way  of  limiting  the  increase  in  carriers  during  this  stage  is  needed.  Ac- 
cording to  Blakemore  C83,  the  band-to-band  Auger  recombination  is  a process 
complementary  to  impact  ionization.  An  expression  of  the  band-to-band 
recombination  C7]  and  £83  is  given  by: 

(122) 

where  B is  a constant.  In  a later  section,  the  computer  simulation  results 
for  different  values  of  B are  presented. 

Near  steady  state  stage 

This  stage  is  described  by  equations  (112)  and  (114),  and  is  shown  in 
Fig.  67.  In  order  to  satisfy  equation  (113)  in  the  case  of  high  current 
density,  the  equilibrium  lifetime  "tq"  must  be  shorter.  Therefore,  the 
equilibrium  lifetime  depends  very  much  on  the  level  of  injection.  An  inves- 
tigation of  the  solution  with  different  equilibrium  lifetimes  is  also  found 
in  a I ater  section. 

The  use  of  transient  I ifetime  and  band-to-band  recombination  in  the  case  of 
an  extreme  current  step 

In  this  section,  the  externally  applied  current  step  is  the  same  as  the 

L.  ? 

one  described  before,  with  10  Amp/ cm  amplitude  and  2 nanoseconds  rise 
time.  The  transient  lifetime  adjustment  described  previously  is  first  used 
without  any  band-to-band  recombination  during  the  impact  ionization  stage. 
The  voltage  vs.  time  transient  response  is  shown  in  Fig.  71  as  the  curve  in- 
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- dicated  by  a zero  band-to-band  recombination  constant  (B  = 0).  Fig.  71  also 

shows  the  transient  response  using  both  the  transient  lifetime  adjustment 
and  band-to-band  recombination  with  different  values  of  the  constant  B.  The 
electron  and  hole  currents,  the  displacement  current,  the  carrier  densities, 
and  the  electric  field  as  a function  of  distance  are  shown  in  Figs.  72  to  75 

for  the  first  part  of  the  transient  response  which  is  common  to  all  curves 

/ 

of  Fig.  71.  At  this  stage  t^e  total  current  in  the  N”-region  is  mainly  car- 
ried by  the  displacement  current  component  as  shown  in  Figs.  72  and  73. 
From  Fig.  74  it  is  clear  that  the  use  of  the  transient  lifetime  does  limit 
the  excess  diffusion  of  carriers  near  the  heavily  doped  regions.  Figures  76 
to  79  demonstrate  the  variations  of  the  physical  quantities  with  no  band- 

I 

to-band  recombination.  As  explained  previously,  the  particle  currents,  as 
well  as  the  carrier  level  increases  to  a very  high  value,  which  cause  the 
displacement  current  to  switch  to  a high  negative  value.  This  results  in  a 
fast  decrease  in  the  electric  field  as  shown  in  Fig.  79  for  the  time  instant 
t^,  and  consequently  the  output  voltage  continues  to  drop  to  a very  low 
value  as  shown  in  Fig.  71.  Figures  80  to  83  demonstrate  the  spatial  varia- 
tion of  the  physical  parameters  using  band-to-band  recombination  with  a 
value  of  5 X 10^^  for  the  constant  B.  At  time  instant  t^g,  the  impact  ioni- 
zation is  talcing  place,  the  electric  field  and  the  carrier  levels  are  high, 
the  sum  of  the  electron  and  hole  currents  exceeds  the  total  applied  current. 
Therefore,  the  displacement  current  switches  to  a negative  value.  Combining 
that  with  the  band-to-band  recombination,  the  net  effect  is  a decrease  in 
the  carrier  levels  and  the  electric  field.  Consequently,  the  electron  and 
hole  curr3nts  decrease,  while  the  displacement  current  switches  back  to  a 
positive  value  as  shown  for  the  curves  and  time  instant  t^^.  At  the  time 
instant  t.2  the  device  reaches  steady  state,  the  displacement  current  goes 


Figure  7i.  V(t)  for  an  extremely  high  current  step  using  the  transient  lifetime  and 
band" to- band  recombination 
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Figure  73.  Displacement  current  density  distribution  for  time  instants 

t,  to  t_  as  shown  In  Figure  71 
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Figure  7^.  Electron  and  hole  density  distributions  for  the  instants 
t * 0 to  tg  as  shown  in  Figure  71 
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Figure  80,  Electron  and  hole  current  density  distributions  for  t 
Instants  t,.  to  t,  as  shown  in  Figure  71 
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I to  zero,  and  the  sum  of  the  electron  and  hole  currents  is  equal  to  the  total 

i 

( current. 

I In  conclusion,  the  use  of  transient  lifetimes  and  band-to-band  recombi- 

1 

, nation,  indeed  predict  a high  voltage  across  P*-N  -N*  structure  for  extreme- 

ly high  current  densities.  The  steady  state  high  voltage  depends  on  the 
strength  of  the  band-to-band  recombination.  The  larger  the  value  of  the 
constant  B is,  the  lower  the  carrier  level.  Therefore,  the  electric  field 
needed  to  support  the  current  is  higher.  The  initial  decrease  in  the  output 
voltage  in  Fig.  71,  followed  by  an  increase  until  steady  state,  is  attribut- 
ed to  the  instantaneous  response  of  the  carriers  to  the  electric  field  as 
discussed  for  the  reverse  biased  case  in  Chapter  VI. 

I 

External  excitation  in  the  form  of  high  current  step 

In  this  section,  the  forward  current  step  has  a 100  amp/cm^  ampi itude 

and  1.25  nanoseconds.  The  effect  of  varying  the  equilibrium  lifetime  on  the 

near  steady  state  stage,  which  is  previously  discussed,  is  presented  here, 

along  with  the  effect  of  using  the  transient  lifetime.  Figure  84  shows  the 

output  voltage  as  a function  of  time  for  three  different  cases.  The  first 

case,  which  is  indicated  by  the  triangular  shaped  marks,  is  the  transient 

response  without  using  the  transient  lifetime,  and  with  an  equilibrium  life- 
— *8 

time  in  the  N -region  of  4 x 10  seconds.  The  electron  and  hole  current 
densities,  the  displacement  current  density,  the  carrier  density,  and  the 
electric  field  as  functions  of  distance  are  shown  in  Figs.  85  to  88  for  the 
time  instants  t » 0,  t^,  ...,  and  t^.  It  should  be  noticed  that  as  steady 
state  is  approached,  and  the  displacement  current  drops  to  a very  small 
value,  electron  and  hole  diffusion  continues  in  order  to  decrease  the  spa- 
tial derivatives  of  JN  and  JP  to  satisfy  equation  (113).  The  time  instant 
t^  of  Fig.  85  shows  that  at  the  ohmic  contact  of  the  P^-region,  30X  of  the 
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Figure  O^i,  V(t)  for  a high  current  step  with  different  lifetimes 
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total  current  is  carried  by  electrons^  while  the  other  70Z  is  carried  by 
holes.  A similar  situation  on  a lower  scale  appears  at  the  ohmic  contact  of 
the  N^-region.  Whereas,  it  is  known  that  the  total  current  at  the  ohmic 
contact  of  the  P^-region  is  totally  carried  by  holes,  and  at  that  of  the 
N*-region  is  totally  carried  by  electrons.  The  second  case  which  is  indi- 
cated by  circular  shaped  marks  in  Fig.  84  illustrates  the  transient  response 
using  the  transient  lifetime  with  the  same  value  of  equilibrium  lifetime  as 
before.  The  spatial  distribution  of  the  important  parameters  for  the  time 
instant  tg,tg,...,  and  t^2  shown  in  Figs.  89  to  92.  As  shown  in  Fig. 
84,  the  use  of  the  transient  lifetime  causes  an  increase  in  the  peak  voltage 
due  to  the  fact  that  it  limits  the  excess  majority  carriers  to  that  of  the 
minority  carriers  in  the  heavily  doped  region,  near  t = 0,  as  explained  pre- 
viously. Due  to  the  equal  equilibrium  lifetimes  of  these  two  cases,  the 
near  steady  state  behavior  is  identical,  as  seen  by  the  equality  of  the  out- 
put voltage  in  Fig.  84.  This  is  also  shown  by  the  carrier  diffusions  con- 
tinuing to  decrease  the  spatial  derivatives  of  JN  and  JP,  as  shown  for  the 
time  instant  t^2  Figs.  89  and  91.  The  third  case  which  is  indicated  by 
the  symbol  "x"  in  Fig.  84  demonstrates  the  transient  lifetime  and  with  an 
equilibrium  lifetime  of  8 x 10  seconds.  This  lifetime  is  50  times  short- 
er than  that  of  the  previous  two  cases.  The  net  effect  of  the  shorter  life- 
time is  an  increase  in  the  steady  state  voltage,  as  well  as  in  the  peak  vol- 


tage. 

Figures  93  to  95  illustrate  the  parameters  of  interest  as  a function  of 
distance  at  steady  state.  The  shorter  equilibrium  lifetime  limits  the  car- 
rier diffusion  and  satisfies  the  steady  state  condition  of  equation  (113) 
which  is  repeated  as  follows: 
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As  a result  the  total  current  in  the  P^-region  is  carried  by  holes  and  by 
electrons  in  the  N*-region  as  shown  in  Fig.  95.  Also  the  lower  carrier  lev- 
els causes  a higher  electric  fields  which  in  turn  results  in  a higher  steady 
state  voltage. 

In  conclusion,  the  results  of  this  chapter  show  that  the  lifetime  of 
the  ShockI ey-Read-Hal I recombination  expression  depends  on  the  level  of  in- 
jection in  the  case  where  the  junction  is  forward  biased.  As  the  injection 
level  increases,  the  lifetime  becomes  shorter  in  order  for  the  device  to 
reach  a steady  state  where  the  total  current  is  carried  by  majority  carriers 
near  the  ohmic  contacts;  and  as  the  I ifetime  becomes  shorter,  the  steady 
state  voltage,  as  well  as  the  peak  voltage  of  the  transient  response,  in- 
creases. 

In  case  of  extreme  levels  of  injection  where  the  device  is  too  far  from 
equilibrium,  the  use  of  a transient  lifetime  is  necessary  in  order  to  limit 
the  excess  diffusion.  If  the  level  of  injection  is  high  enough  to  cause 
avalanche  multiplication,  it  becomes  necessary  to  use  band-to-band  Auger 
recombination  to  limit  the  plasma  in  the  avalanching  region. 


I 


i. 
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CHAPTER  VI 

EXPERIMENTAL  RESULTS 

This  chapter  outlines  the  diode  fabrication  procedures  as  well  as  the 
methods  employed  to  determine  the  important  diode  parameters.  These  parame- 
ters are  the  diode  area  (A)^  the  N”-region  doping  density  <Np),  and  its 
width  (H) . These  diodes  were  tested  in  the  forward  direction  with  an  exter- 
nal excitation  in  the  form  of  a current  pulse.  A comparison  between  the  ex- 
perimental results  and  the  computer  simulation  is  also  included. 

As  a point  of  comparison,  N|^  was  first  determined  by  calculations  based 
on  the  depletion  C-V  characteristics  of  MOS  capacitors,  which  were  fabricat- 
ed from  the  same  starting  wafer.  Then  N^  was  determined  from  electrical 
measurements  on  fabricated  diodes.  The  diodes  are  diffused  and  were  fabri- 
cated with  very  short  drive  times  (two  hours  at  1100°C)  for  the  boron  (P*) 
diffusion.  This  procedure  results  in  diodes  with  junctions  that  behave  as 
nearly  abrupt  junctions.  This  was  verified  when  the  slope  of  the  reverse 
C-V  characteristics  (on  log-log  paper)  approached  -0.5.  Once  the  abrupt 
junction  was  confirmed,  the  theory  available  allowed  calculation  of  N^  from 
the  slope  of  versus  V characteristics. 

The  width  of  the  n”  region  of  the  diode. was  measured  directly  by  angle 
lapping  and  then  staining  with  an  acid  solution.  The  wafer  was  scribed  and 
broken  to  get  a chip  of  the  proper  size  for  mounting  on  the  5°  chuck  for 
lapping.  Once  stained,  the  calibrated  eye  piece  in  the  microscope  allowed 
measurement  of  the  width  W.  To  get  the  width  U from  electrical  measure- 
ments, a doping  profile  is  required.  This  can  determine  the  width  of  the 
lightly  doped  region  if  punch  through  precedes  impact  ionization  in  the  re- 


versed  biased  direction,  as  shown  for  the  first  device  considered  in  Chapter 
IV.  If  avalanche  occurs  first,  then  the  N”-region  is  not  completely  deplet- 
ed, and  W cannot  be  determined  from  the  C-V  measurements. 

The  specification  for  the  starting  wafers  of  single  crystal  silicon  are 
listed  in  Table  4.  All  the  wafers  were  <111>  oriented  epitaxial  N on  N*. 
The  resistivity  and  width  values  listed  were  from  data  provided  by  the 
manufacturer. 

Table  4.  Specifications  of  the  Silicon  Wafers 


Wafer 

Width 

Resistivity 

Doping 

(#} 

(um) 

(ohm-cm) 

(cm 

1 

12.10 

8.13 

6.0  X lo"*^ 

2 

52.90 

28.60 

1.8  X lO**^ 

3 

21.34 

2.48 

2.0  X lo”*^ 

4 

54.80 

29.30 

1.7  X 10^*^ 

The  diode  fabrication  procedures  as  wel I as  the  MOS-C  procedure  are 
outi ined  bel ow: 

Diode  Fabrication  Procedure 

1.  The  wafer  was  ultradeaned  by  running  through  a series  of  solvents  and 
acids. 

2.  A 2000X  thick  oxide  was  grown  on  the  wafer  in  a dry  oxygen  atmosphere  at 
1100°C  for  4 hours  and  30  minutes. 
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3.  The  wafer  was  scribed  in  hal  f to  be  processed;  one  as  diodes  and  the  oth- 


er as  NOS-capacitors. 


4.  By  using  photoresist  and  masking  procedures^  the  diffusion  windows  were 
opened  in  the  oxide. 


5.  A boron  deposition  of  20  minutes  at  1100°C  was  followed  by  a 2-hour  drive 
in  a dry  oxygen  atmosphere  at  1100°C. 

6.  The  oxide  stabilization  step  was  performed  in  the  phosphorus  deposition 
furnace  for  10  minutes. 


7.  A photoresist  masking  procedure  was  used  to  open  contact  windows  in  the 
oxide. 


8.  This  was  followed  by  an  aluminum  evaporation  step  at  10*°  torr^  deposit- 


ing a coat  of  2000{(  to  3000i{  of  aluminum. 


9.  The  photoresist  masking  procedures  were  used  to  remove  the  excess  alumi- 
num and  to  form  the  desired  aluminum  pattern  on  the  wafer. 

10.  An  aluminum  evaporation  step  was  performed  on  the  back  of  the  wafer, 
after  stripping  the  oxide  while  protecting  the  wafer  front. 


MOS-C  Fabrication  Procedure 


1.  A phosphorous  gettering  step  was  performed  in  the  phosphorus  deposition 
furnace  for  20  minutes  after  the  oxidation  step. 


2.  The  aluminum  evaporation  step  was  performed  at  10  torr,  depositing  a 
coat  of  2000S  to  3000ii  of  aluminum. 
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2.  The  wafer  was  annealed  for  5 minutes  at  500°C  in  a nitrogen  atmosphere. 


4.  The  last  two  steps  of  the  diode  fabrication  were  then  repeated  to  com- 
plete the  MOS-C  fabrication. 


N,  from  MOS  Measurements 


To  insure  that  the  starting  material  was  correct  as  labeled^  MOS  capa- 
citors were  fabricated  so  that  the  epitaxial  layer  doping  density  could  be 
established  before  the  other  methods  were  attempted.  Fig.  96  illustrates 
the  MOS  structure. 

For  the  MOS  wafers,  the  doping  density  was  found  by  obtaining  the  accu- 
mulation capacitance  and  the  minimum  depletion  capacitance  from  the  C-V 
characteristic  similar  to  Fig.  97.  The  theory  for  this  is  presented  below. 


. Cf,  = capacitance  of  the  SiO^ 

r - “ 

''equiv.  “ C.+C-  C-  = capacitance  of  the  semiconductor 

S U d 


For  the  geometry  of  Fig.  96,  we  have  the  series  combination  of  the  two  capa- 
citances. 


equiv. 


1 + Cq/Cj  *^0“ 

^ * XqCqKjA  '' 


•^S^o  * *^0^ 


^oVs  ' ^eq'^S^O  ‘^eq'^o“ 
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^ Ws  - ^ ^ fio-  - 

'eq'^0  >^0  [‘^eq  J 


and  therefore 


[c^-’l 

L imn  J 


(123) 


Cq  is  the  maximum  capacitance  of  accumulation  where  the  voltage  is  positive. 
In  accumulation  = Cg.  Having  Cg^  the  width  of  the  oxide  layer  is  cal cu- 
I ated. 


(124) 


This  value  is  a much  more  accurate  method  for  determining  Xg  than  might  be 
estimated  from  the  color  code  for  silicon  dioxide  on  silicon. 

Then^  U was  calculated,  where: 

indX 

K = 1.38  X 10~^^  Joule/°K 


T = 300°K 

Kg  = 11.8 

Eg  = 8.854  X 10  ^ Fd/cm 
q = 1.602  X 10“^’  coul 
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n.  = 1.18  X 10*'°  cm"^ 


Kg  = 3.8 


A = 12.6  X 10  ^ cm^ 


Having  W ^ was  obtained  from  a plot  made  from  the  following  formula  C6] 
(fl3X  1/ 

relating  N.  and  W_,  . 

D mau 


1 

max 

^4<T  Kg  Cq 

Nk 

'"TT 

"i 

Table  5 summarizes  the 

data  found  for  the 

MOS  wafers. 

Table 

5.  MOS 

Anal ysis  of 

Njj  for  Starting  Materia' 

Wafer 

C 

C . 

Xn 

N,  (MOS) 

(rated) 

max 

min 

0 

D 

D 

n 

(pf) 

(pf) 

(A) 

(cm)  ^ 

(cm)  ^ 

1 

12.4 

1.4 

3420 

3.60x10*'^ 

6.0x10*'^ 

2 

11.2 

4.2 

3790 

1.65x10^^ 

1.8x10‘'^ 

3 

10.9 

6.8 

3890 

15 

1.75x10'^ 

15 

2.0x10 

4 

16.2 

5.6 

2620 

2.60x10^^ 

1.7x10*'^ 

(125) 


^ Diode  £~V  Measurements 

For  the  diode  wafers^  the  reverse  C-V  characteristics  and  the  log  C vs. 
log  V of  these  characteristics  were  plotted  as  shown  in  Figs.  98  and  99  for 
wafer  #3.  This  was  done  to  see  if  the  junction  can  be  approximated  by  an 


A ^ 


A0-A0a3  539  PURDUE  UNIV  LAFAYETTE  IND  SCHOOL  OF  ELECTRICAL  ENGI— ETC  F/6  9/1 

THE  ANOMALOUS  VOLTAGE  RESPONSE  OF  THE  P(+)N-N(+)  DEVICE  AND  ITS~ETC(U) 
DEC  77  6 W NEUDECKr  L R RAZOUK  DAAe29-77-6-0191 

UNCLASSIFIED  TR-EE-77-40  ARO-15322. 1-A-EL  NL 


3of3 

AD 

a 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

■ 

■ 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

END 

OATC 

riLMO 

6-78 

DOC 

Capacitance  In  PF 


abrupt  junction.  An  abrupt  junction  Is  known  to  have  a log  C vs.  log  V plot 
with  a -0.5  slope.  If  the  junction  Is  linearly  graded,  this  slope  becomes 
-0.333  C63,  These  slopes  were  found  to  be  -0.29,  -0.5,  -0.51,  and  -0.42, 
for  wafers  #1,  #2,  #3,  and  #4  respectively.  For  the  abrupt  junction  diodes, 
the  reverse  biased  C vs*  voltage  can  be  used  to  determine  Nj^.  In  this  case, 
the  following  equations  are  applicable: 


(126) 


q = 1.602  X 10’^’ 

for  the  semiconductor 

A = area  of  the  region 


= built  In  potential,  typically  less  than  1 volt. 


1 .2<“ba*> 


F 


(127) 


(128) 


Equation  (128)  Indicates  the  slope  of  the  1/C^  vs.  (Vj^^V)  plot  and  Its 
relationship  with  N|^.  For  a uniformly  doped  N~  region,  this  slope  would  be 
constant.  For  a non-uniformly  doped  N~  region,  the  slope  would  change  and 


2 

the  plot  is  not  a straight  line.  Fig.  100  illustrates  the  1/C  vs.  V plot 
for  wafer  #3.  The  plot  is  linear.  Therefore,  the  background  doping  is  uni> 
forn  for  this  wafer. 

The  doping  profile  method  consists  of  plotting  the  1/C^  vs.  V charac- 
teristics of  tha  wafer  all  the  way  out  to  the  breakdown  voltage.  N|^(W)  is 
found  in  the  same  manner  as  (diode).  In  this  case,  however  N^^CW)  is 
found  from  the  slopes  of  adjacent  points  on  the  1/C^  vs.  V characteristics. 
The  width  W at  each  point  is  found  from  the  abrupt  junction  capacitance 
equation  as  follows: 


Table  6 summarizes  the  results  for  N^.  In  obtaining  these  results  the 
abrupt  junction  approximation  was  used  for  all  wafers. 


Table  6.  Results  for  Nj^  Using  Different  Methods 


Wafer 

(diode) 

N|^  (rated) 

N|^  (NOS) 

It 

(cm)”^ 

(cm)”^ 

(cm)"*^ 

1 

6.0  X lO**^ 

3.60  X lO**^ 

4 

2.67  X 10^* 

1.8  X lo"'^ 

1.65  X 10^^ 

5 

2.64  X 10**^ 

2.0  X 10^^ 

1.75  X 10**^ 

6 

3.46  X 10^^ 

1.7  X 10^^ 

2.60  X 10^^ 

Note  that  the  Nj^  (diode)  was  higher  than  the  rated  value  and  that  NOS 
data  gave  a value  less  than  the  diode  data.  This  seems  reasonable  since 
doping  with  impurities  at  high  temperatures  during  the  diffusion  for  the 
diodes  may  add  unwanted  impurities. 


urenents.  If  the  breakdown  Mechanism  is  pure  punch-through^  then  U could  be 
found  from  the  Np(W)  profiling  technique.  The  more  usual  case^  however,  is 
pure  avalanching. 

A destructive  test  to  obtain  U is  illustrated  in  Fig.  101.  The  same 
angle  lap  and  stain  procedure  is  used  to  get  the  junction  depth  (of  the 
region)  and  the  N**  region  W. 

The  major  effort  for  measuring  these  parameters  consists  of  lapping  and 
staining.  Chips  were  cut  to  the  proper  size  for  mounting  on  the  5°  chuck, 
and  the  chips  were  angle-lapped  until  the  proper  depths  had  been  reached. 
Then,  the  lapped  regions  were  stained  with  the  hydrofluoric  acid  solution. 
It  is  often  necessary  to  repeat  the  lapping  and  staining  several  times  in 
order  to  clearly  discern  the  interfaces  between  the  regions  of  the  wafer. 
Since  a 5°  chuck  was  used,  the  geometry  of  Fig.  101(b)  is  applicable.  The 
microscope  can  be  calibrated  such  that  a calibrated  vernier  is  read  through 
the  microscope  and  knowing  the  lens  magnification  the  microscope  measured 
depth  is  easily  determined. 

Actual  depth  > (sin  5*^) (microscope  meas.  depth)  (130) 

Then,  by  equation  (130),  the  depth  of  either  the  P'*’  or  n“  region  is  calcu- 
lated. It  is  recommended  that  for  a shallow  region,  a 1°  block  be  used 
while  for  a large  W a 10°  block  be  used  to  ease  the  angle  lapping  and  stain 
problems. 

Table  7 shows  the  results  of  determining  the  P^  region  and  n'  region 
widths  by  the  angle  lapping  and  stain  procedures. 


5*  Angle  Block 


Figure  101.  Geometry  of  Angle  Lap  and  Stain 


Samples  of  the  fabricated  diodes  from  wafers  M2,  MZ,  and  #4  were  pack- 
aged and  tested  with  large  forward  current  pulses.  The  output  voltage  of 
the  Velonex  pulse  generator  was  applied  to  the  series  combination  of  diode 
and  a 225  ohm  resistor.  This  results  in  a current  pulse  through  the  diode. 
The  Velonex  pulse  generator  has  a maximum  output  voltage  of  2.5  kilo  volts, 
with  a 30  to  50  nanosecond  rise  time.  This  rise  time  is  a function  of  the 
load  impedance  which  results  from  the  resistor  and  the  diode  series  combina- 
tion. Current  pulses  of  0.48,  0.80,  1.60,  and  3.20  Amperes  amplitude,  50 
nanoseconds  rise  time,  and  300  nanoseconds  pulse  width,  were  used  to  obtain 
the  diode  transient  response.  The  current  through  the  diode,  and  the  vol- 
tage across  it  were  recorded  as  functions  of  time  from  a dual  trace  oscillo- 
scope. Table  8 summarizes  the  results,  where  T^^,  Tj^p,  Vj^p,  and 

are  shown  in  Fig.  102  along  with  the  general  shape  of  the  transient 
response.  The  quantities  of  Table  8 are  defined  as  follows: 

T is  the  rise  time  of  the  applied  current  pulse 
rs 


i 


J 


I 


Figure  102,  External  applied  current  and  output  voltage 
for  experimental  data 
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t : 

I Tjjp  is  the  tine  corresponding  to  the  peak  voltage  across  the  diode 

Vjjp  is  the  peak  voltage  across  the  diode 
Vdss  ■’s  steady  state  voltage  across  the  diode 

j 

2 

I are  the  amplitude  of  the  current  pulse  in  ampers  and  amp/cm 
D V j 

' 1 

Table  8.  Results  of  pulsing  the  diodes  in  the  forward  direction  | 


Wafer  # 

Tps 

^DP 

'^OP 

''oss 

(n.sec.) 

(n.sec.) 

(vol ts) 

(volts) 

(amp) 

(amp  cffl^) 

2 

50 

50 

11.2 

4 ' 

0.48 

3.81  X 10^ 

2 

50 

50 

13.4 

5 

0.80 

6.35  X 10‘ 

2 

50 

40 

23.5 

7 

1.60 

1.27  X 10^ 

2 

50 

40 

36.0 

10 

3.20 

2.54  X 10^ 

3 

50 

40 

3.4 

2.0 

0.48 

3.81  X 10^ 

3 

50 

40 

4.6 

2.2 

0.80 

6.35  X 10^ 

3 

50 

40 

7.6 

3.2 

1.60 

1.27  X 10^ 

3 

50 

40 

15.0 

5.0 

3.2 

2.54  X 10^ 

4 

50 

30 

5.5 

2.6 

0.48 

3.81  X 10^ 

4 

50 

30 

7.3 

3.6 

0.80 

6.35  X 10^ 

4 

50 

30 

10.8 

4.6 

1.60 

1.27  X 10^ 

4 

50 

35 

18.5 

7.0 

3.10 

2.46  X 10^ 

The  diodes  were  also  pulsed  in  the  reverse  biased  direction.  The  pulsed  I-V 
characteristic  is  shown  in  Fig.  103  for  wafer  #3. 


Comparison  between  Theory  and  Experiment 

Diodes  of  wafer  #3  were  chosen  for  comparison  with  the  computer  Simula-* 
tion  results  under  reverse  and  forward  pulsed  conditions.  A summary  of  the 


J 
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measured  physical  parameters  for  diodes  of  wafer  #3  are  as  follows: 

Type  of  junction  = abrupt 

Doping  density  = 2 x 10^^/cm^  rated  value 

= 1.75  X 10^^/cm^  MOS-C  measurements 
15  3 

= 2.64  X 19  /cm  diode  C-V  measurements 

Depth  of  junction  =1.4  microns 

N”-region  thickness  = 17  17  microns 

-4  2 

Diode  area  = 12.6  x 10  cm 

Reverse  Pul sed  Case 

The  doping  density  of  the  N”-region  was  first  assumed  to  be  2 x 
10  /cm  according  to  the  rated  value,  because  it  falls  between  the  two 
values  obtained  from  the  HOS-C  and  the  diode  C-V  measurements. 

A transient  behavior  solution  similar  to  the  one  discussed  in  Chapter 
IV  was  obtained.  The  external  excitation  was  in  the  form  of  a current  step 
with  a rise  time  of  0.1  nanoseconds,  and  amplitude  of  100  Amp/cm^,  or  126 
milli-Amp.  The  computer  simulation  resulted  in  a steady  state  output  vol- 
tage of  185  volts.  This  voltage  is  very  high  compared  to  the  experimental 
data  of  the  pulsed  I-V  characteristic  of  Fig.  103.  This  indicated  that  the 
actual  doping  density  of  the  N -region  is  higher  than  the  rated  value. 
Therefore,  the  value  of  obtained  from  the  diode  C-V  measurements  of  2.64 
X 10^^/cm^  was  assumed.  The  transient  behavior  under  reverse  pulsed  condi- 
tions for  this  case  is  included  in  Chapter  III  under  the  same  external  exci- 
tation previously  discussed  of  126  milli-Amp  amplitude.  The  result  of  the 
output  voltage  as  a function  of  time  is  shown  in  Fig.  39  which  shows  that 
'the  steady  state  voltage  is  150  volts.  This  value  is  in  excellent  agreement 
with  the  experimental  data  of  Fig.  103.  Therefore,  the  actual  doping  densi- 
ty of  the  N“-region  is  that  obtained  from  the  diode  C-V  measurements. 
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According  to  Sze  C63^  an  approximate  expression  of  the  breakdown  vol- 
tage for  a one  sided  abrupt  junction  is  given  by  the  following  equation: 

Vg  = 60(Eg/1.1)^^^  volts  (131) 

where  Eg  is  the  band  gap  in  eV(=1.12  eV  for  Si).  Equation  (131)  results  in 
a breakdown  voltage  of  167.3  volts  which  is  higher  than  the  actual  value  of 
150  volts.  It  should  also  be  noted  that  the  approximate  equation  of  (131) 
predicts  a constant  breakdown  voltage  for  all  current  levels^  while  the 
steady  state  condition  of  the  computer  simulation  depends  on  the  current 
level  according  to  equations  (106)  and  (107)  of  Chapter  IV. 

Chapter  III  presented  the  thermal  equilibrium  condition  for  this  device 
assuming  a total  device  length  of  20  microns.  The  importance  of  using  non- 
uniform  spatial  steps  is  also  discussed  in  Chapter  III. 

Forward  Pul sed  Case 

The  first  part  of  the  measured  transient  response  under  forward  pulsed 
conditions^  which  is  described  in  Table  8^  is  shown  in  Fig.  104  for  diodes 
of  wafer  UZ,  with  a pulse  amplitude  of  3.2  Amp.  (or  2.54  x 10^  Amp./cm^)/ 
and  a rise  time  of  50  nanoseconds.  One  might  be  tempted  to  try  to  simulate 
this  behavior  with  a current  step  of  long  rise  time  as  used  experimentally. 
This  proved  to  be  impractical  for  several  reasons.  First,  it  should  be  not- 
ed that  a computer  simulation  does  not  contain  acceleration  terms,  the  car- 
rier response  to  the  electric  field  is  essentially  instantaneous.  This 
results  in  a shorter  response  time  when  compared  with  experimental  data. 
For  example,  if  the  measured  peak  voltage  occurs  at  time  Tj^p,  the  computer 
simulation  results  in  a peak  voltage  at  a much  shorter  time.  Therefore,  the 
computer  simulation  will  predict  and  explain  the  device  behavior  with  a dif- 
ferent time  scale.  As  a result  the  use  of  long  rise  times,  as  used  experi- 
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Figure  10*I.  V(t)  obtained  experimentally  for  diode  #3 


■entaMy,  will  not  yield  an  exact  response.  It  is  also  inportant  to  realize 
that  the  tiae  step  required  to  achieve  an  accurate  converged  solution 
depends  very  much  on  the  size  of  the  saallest  spatial  step^  as  well  as  the 
shortest  lifetime  used.  This  tiae  step  limitation  was  analyzed  in  detail 
for  the  case  of  explicit  method  Z21,  while  due  to  the  complicated  nature  of 
the  implicit  method^  a relation  between  the  required  time  step  and  the  dif~ 
ferent  parameters  used  is  hard  to  find.  Therefore^  the  use  of  a non-uniform 
spatial  mesh  with  small  steps  near  the  P^-N  junction  limits  the  time  step 
size.  AlsO/  the  use  of  a non-equilibrium  transient  lifetime  adds  another 
restriction  to  the  time  step  size,  due  to  the  presence  of  very  short  life- 
times in  the  beginning  of  the  transient  solution.  Therefore,  in  order  to 
reduce  the  computational  time  to  a practical  level  a short  rise  time  is  used 
in  general.  This  rise  time  is  further  reduced  if  a transient  lifetime  is 
used. 

As  discussed  in  Chapter  V,  the  Schockl ey-Read-Hat I lifetime  depends  on 
the  level  of  injection.  A shorter  equilibrium  lifetime  results  in  a high 

m 

steady  state  voltage,  while  the  use  of  transient  lifetimes  results  in  a 
higher  peak  voltage. 

t 

Figures  105  and  106  show  the  transient  response  of  the  computer  simula- 
tion for  the  experimental  device  with  the  external  excitation  in  the  form  of 

3 2 

a current  step  with  an  amplitude  of  2.54  x 10  Amp/ cm  . Figure  105  shows 
the  result  with  a rise  time  of  1 nanosecond,  and  an  equilibrium  lifetime  of 
3.1  X 10~^^  seconds.  This  results  in  a peak  voltage  of  11  volts  and  a 
steady  state  voltage  of  7.6  volts.  It  is  clear  that  the  steady  state  vol- 
tage is  higher  than  the  experimental  value  which  implies  that  a longer  life- 
time is  required,  while  the  peak  voltage  which  is  lower  than  the  et  'erimen- 
tal  value  implies  the  need  of  using  the  transient  lifetime.  Therefore,  the 


volts 
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simulation  was  repeated  using  the  transient  lifetime  adjustment  and  an 
equilibrium  lifetime  of  4 x seconds;  also  the  rise  time  was  decreased 
to  1 picosecond.  Figure  106  shows  the  result  of  this  transient  behavior. 
The  peak  voltage  has  a value  of  17.4  volts  compared  to  an  experiaental  value 
of  15.7  volts,  while  the  steady  state  voltage  is  6.8  volts  compared  to  an 
experimental  value  of  5 volts  which  implies  that  the  actual  equilibrium 
lifetime  is  longer  than  the  used  value  of  4 x 10  seconds.  A longer  life- 
time will  result  in  a lower  steady  state  voltage  as  well  as  a lower  peak 
value.  In  other  words  approach  the  experinental  results. 

Therefore,  the  computer  simulation  can  be  used  to  obtain  a value  for 
the  equilibrium  lifetime  in  the  forward  biased  direction  by  comparing  the 
computer  simulation  results  with  the  experimental  results. 
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CHAPTER  VII 

SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY 
Suwary  and  conct usions; 

A suMHiary  of  the  research  presented  in  this  report  and  the  conclusions 
which  aay  be  drawn  fro«  this  research  are  given  in  the  present  chapter.  On 
the  basis  of  the  results  of  this  work,  some  recommendations  are  then  made 
for  further  study. 

Numerical  methods  were  used  to  solve  the  basic  one-dimensional  two  car- 
rier transport  equations  coupled  with  Maxwell's  and  Poissons's  equations. 
These  solutions  were  used  to  analyze  and  to  explain  the  transient  behavior 
of  the  P*-N-N^  structure,  in  case  of  extremely  high  injection  levels,  for 
both  the  reverse  and  forward  directions. 

Special  effort  was  directed  toward  explaining  the  anomalous  voltage 

that  appears  across  the  diode  terminals  when  pulsed  in  the  forward  direction 

5 2 6 2 

with  an  intense  current  density  of  10  amp/ cm  to  10  amp/ cm  . The  problem 
formulations  described  are  of  a very  general  character  with  none  of  the  con- 
ventional assumptions  and  restrictions  introduced.  There  is  freedom  in  the 
choice  of  the  doping  profile,  generation-recombination  law,  mobility  depen- 
dencies, and  injection  level.  The  boundary  conditions  are  applied  at  the 
external  ohmic  contacts.  For  a specified  arbitrary  current  excitation,  the 
solution  yields  the  terminal  voltage  and  all  the  quantities  of  interest 
throughout  the  entire  device  as  functions  of  position  and  time. 

A comparison  between  the  explicit  and  the  implicit  numerical  method 
shows  that  although  the  implicit  method  is  complicated  and  hence  is  more  ex- 
pensive for  each  time  step  At,  it  is  more  stable  numerically  and  the  time 
step  At  can  usually  be  made  very  large  when  the  solution  approaches  steady 
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state  conditions.  Since  the  implicit  technique  requires  a good  starting 
point,  the  explicit  technique  is  used  to  obtain  a first  estinate  for  the  i»- 
pMcit  technique.  An  accurate  Solution  of  the  theraal  equilibrium  condition 
problem  is  also  essential  in  obtaining  a convergent  solution  of  the  tran- 
sient behavior  at  extreme  injection  levels. 

A numerical  method  was  used  to  obtain  an  accurate  distribution  for  the 
electric  field,  the  hole  density,  and  the  electron  density  at  thermal 
equilibrium.  In  order  to  numerically  obtain  a zero  value  for  both  the  elec- 
tron and  the  hole  current  densities,  different  space  charge  points  and  field 
points  were  used.  The  thermal  equi ibrium  results  show  that  non-uniform  spa- 
tial steps  are  required  in  the  case  of  large  potential  variations  across  the 
junction. 

It  was  determined  that  devices  with  a low  doping  density  and  a short 
N~-region  deplete  at  a relatively  fast  rate  until  punch  through  under 
reverse-pulsed  conditions.  The  electric  field  increases  uniformly  in  the 
N”-region,  until  the  onset  of  avalanche  by  impact  ionization.  On  the  other 
hand,  a device  with  high  doping  density  and  long  N -region  does  not  fully 
deplete.  In  this  case,  the  electric  field  increases  while  keeping  its  tri- 
angular shape  in  the  depletion  region  until  it  reaches  avalanche.  In  both 
cases,  the  steady  state  condition  is  satisfied  when  the  rate  of  generated 
carriers  via  avalanche  multiplication  is  balanced  by  the  recombination  and 
depletion  rate. 

For  the  forward  pulsed  case,  the  use  of  the  standard  Shocki ey-Read-Hal I 
recombination  term  is  adequate  in  explaining  the  transient  response  for  re- 
latively low  current  densities,  but  it  fails  to  give  a physical  explanation  1 

for  extreme  current  levels,  where  the  device  is  too  far  from  the  eqil ibrium 
condition.  This  is  mainly  due  to  the  lack  of  sufficient  recombination  dur- 
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ing  the  transient  behavior.  The  transient  problem  can  be  divided  into  three 
stages,  the  excess  diffusion  stage  near  t = 0,  the  impact  ionization  stage, 
and  the  near  steady  state  stage.  In  the  first  stage,  the  use  of  the  tran- 
sient lifetimes,  which  are  adjusted  in  such  a way  as  to  limit  the  increase 
of  majority  carriers  to  that  of  minority  carriers  in  the  two  heavily  doped 
regions,  is  necessary  in  order  to  limit  the  excess  diffusion.  These  tran- 
sient lifetimes  reduce  to  their  equilibrium  value  as  steady  state  is  ap- 
proached. In  the  impact  ionization  stage,  the  band-to-band  Auger  recombina- 
tion is  necessary  to  limit  the  plasma  formation  in  the  avalanched  region. 
The  behavior  of  the  near  steady  state  stage  necessitates  the  use  of  a short 
equilibrium  lifetime  which  depends  on  the  level  of  injection  to  insure  that 
the  total  current  is  carried  b>  tne  majority  carriers  in  the  two  heavily 
doped  regions. 

The  results  show  that  the  use  of  transient  lifetimes  and  band-to-band 
Auger  recombination  indeed  predicts  a high  terminal  voltage  across  the  diode 
when  pulsed  with  intense  current  density  in  the  forward  direction.  A com- 
parison between  experiment  and  theory  shows  that  the  computer  simulation  can 
be  used  to  explain  the  anomalously  large  voltages  present  in  the  experimen- 
tal data. 

The  effects  of  this  anomalous  voltage  response  that  can  occur,  in  the 
forward  and  reverse  pulsed  direction,  on  second  breakdown  are  as  follows: 


1.  In  the  forward  pulsed  case  the  larger  than,  expected  response,  at  extreme 
currents  at  very  short  pulse  times,  indicates  a larger  energy  input, 
hence  a shorter  time  for  second  breakdown  to  occur  than  would  be  expect- 
ed. These  conditions  would  only  occur  for  very  short  intense  current 


pulses  in  the  nanosecond  regime  and  could  predict  an  energy  input  of  1C 
to  150  times  greater.  The  net  result  is  to  effect  (shorten)  the  time  dc- 
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lay  before  second  breakdown  occurred. 

2.  In  the  reverse  pul sed  case  the  diode  that  breaks  down  due  to  punch  | 

through  has  little  if  any  anomalous  voltage  and  therefore  would  act  ex- 
actly  the  same  as  the  normal  long  current  pulse  case  as  far  as  second  I 

breakdown  is  concerned.  I 
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3.  The  reverse  pulsed  case  with  a normal  avalanche  breakdown  will  have  an 
excess  voltage  that  exceeds  the  normal  breakdown  by  up  to  a factor  of  5 
to  8.  Again  this  will  cause  the  energy  input  to  exceed  the  expected 
value  and  initiate  second  breakdown  prematurely. 

Recommendation  for  further  study: 

The  computer  simulation  presented  in  this  report  is  a very  powerful 
tool.  It  can  be  modified  to  study  the  transient  behavior  between  thermal 
equilibrium  and  steady  state  or  between  two  steady  states,  for  any  one- 
dimensional,  two-terminal  semiconductor  structure.  The  simulation  can  also 
be  modified  to  analyze  the  small  signal  behavior  of  microwave  diodes  CIO]. 

This  research  can  be  extended  by  performing  more  theoretical  and  exper- 
imental comparisons  for  different  diode  structures  and  injection  levels  in 
the  forward  direction.  This  may  lead  to  relationships  describing  the  depen- 
dence of  the  equilibrium  lifetime  and  the  diode  parameters,  such  as  the  dop- 
ing density  and  the  length  of  the  N”-region,  and  on  the  level  of  injection, 
on  the  transient  lifetime  and  Auger  recombination  effects. 

Also  the  forward  transient  behavior  can  be  explained  from  a physical 
point  of  view  if  one  can  measure  the  lifetimes  during  the  transient  periods 
near  t = 0. 
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Due  to  the  limitations  of  our  pulse  generator  and  equipment  the  full 
range  of  the  theory  could  not  be  tested.  Such  tests  should  be  conducted. 
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It  would  also  be  of  interest  to  try  these  siaulations  and  experiaents  on 

silicon-on-sapphire  devices  and  aake  direct  coaparisons  of  channeling  and 

; 

second  breakdown  with  the  excess  energy  input. 
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APPENDIX  A 


THE  FORMULATION  FOR  THE  CONSTRUCTION  OF  THE 


MATRIX  ELEMENTS  IN  THE  IMPLICIT  METHOD 


The  terms  in  the  matrix  elements^  of  equations  (44)  to  (46)  can  be  cal 
cut ated  as  fol I ows: 

For  Ed)  > 0 
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36(1) 

afni 


a6(l) 

■jimy 


» - VNEd) 


• N(J)  - 


VN(I) 

mr 


aE(I) 


Ax 


CN(J2)  - N(J)3 


» - VN(1)  + 


1_ 

AX 


VN(I) 

Trrr 


. 1 IVNCDI 

■ lEd)  I 

= - VPd)  • P(J2)  + ll^l  • ^ CP(J2)  - P(J)D 

jlypci)! 

= - VPEd)  • P(J2)  + L . 1-  CP(J2)  - P(J)D 


1 VPCI) 

AX  E^I) 


= - VPCI)  + ^ 


VPd) 

inr 


ALPd)  • 

= ALPEd) 

+ ALPd) 

» ALPd)  • 


|JNd)|  ♦ BETCI)  • 

• |JNd)|  + BETECI) 
aJNd)  , JNd) 

* imr  rjyrriTT 

ajNd)  , JNCi) 

mrr  UMrnr 


|jpd)l 

• |jpd)l 

+ BETd)  • 


ajpd)  . 

aFTTf 


jpd) 

|jp<i)| 
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3G(I) 


JPUT 


= BET(I)  • 


aG(i) 

aN<J2) 


= AU»(I)  • 


3G(1) 

iPdji) 


= BET(I)  • 


3P(J) 

|jpa)| 

. 3JN(1) 

. JN(1) 

3N02) 

wttt 

. 3JP(1) 

. jp(i) 

_ N(J)  • P<J)  - 1 

* T'cT  ■+  Nurv-^in 


3R(J)  _ (P(J)  ♦ 2)P(J)  + 1 


TqC2  ♦ N(J)  + P(J)3^ 


3R(J)  _ (N(J)  ♦ 2)N(J)  ♦ 1 

TqC2  + N<J)  ♦ P<J)D^ 


For  band-to-band  reconbi nation; 


RggCJ)  = B • CP(J)  • N<J)  - ID 


3Rbb(J) 


3P(J) 


= B • N(J) 


3Rbb(J) 


= B • P(J) 


3N(J) 
where^ 

VN(I)  » VNCECDD 


VNE(l) 


3VN(I) 


3ETrr 


f 


208  - 


VP(I)  = VPCE(I)3 


VPE(I) 


aVP(I) 

TeTiT" 


ALP(l)  - aCE(I)D 


ALPE(I)  = 


3ALP(I) 

aE(I) 


BET(I)  = 0CE(I)3 


BETE(I)  = 

BETEa)  jgjjj 
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APPENDIX  B 

THOMAS  ALGORITHM  FOR  TRI DIAGONAL  MATRIX  C9] 


The  equations  are: 

« b,u,  ♦ = d, 


for  1 < i < R 


with  a^  = C|^  = 0 


The  algorithm  is  as  follows: 
First,  compute 


.... 


»i  ' *1  ' "i 


1 


- 210  - 


APPENDIX  C 

ALGORITHM  FOR  TRI-TRIDIAGONAL  MATRIX  Z91 


The  equations  are 


a">u.  . ♦ a<2>v.  . a a«>».  , ♦ * b'^'v,  * b”>y, 

1 1-1  1 1-1  1 1-1  t ■»  11  11 


a<‘>u.  . » a">,.  , ♦ a‘‘>..  , ♦ b<‘'u.  * b">».  * b<".< 
1 1-1  1 1-1  1 1-1  11  11  11 


+ + C.  + c. 


= d'2) 


♦ ^•®^''i-i  + ^•’^^-1  ♦ + b\^\.  * b$’^. 


4.  <7)..  ^ ,(8)„  . ^(9)^ 

♦ c.  ♦ c. 


bP> 


for  1 < 1 < R 


with  a^^"^  = =0  for  1 < m < 9 


The  algorith*  is  as  follows: 
First  coMpute 


(3)  (7) 
®i  *i-1 


o<1)  - h<1>  - _ a<2).(4) 

Bi  - b.  - a.  a^  x^_, 

.<2)  , h<2)  . *<1>.<2)  - a<2).(5)  _ .(3) 

Bi  = b^  - a^  x^_^  a^  x^.^  a. 


'■■1 

^.1 


11 
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(3) 

-(2) 

,<6) 

.(3) 

,<9) 

'i 

‘i-1  ’ 

*1-1 

- *1 

*1-1 

.(4) 

(1) 

.(5) 

,<4) 

,<6) 

,<7> 

"i-1 

*1-1 

“ «1 

*1-1 

(2)  . (5) 

^-1  ®1 

X^5> 

*1-1 

- a^*^ 
®1 

*1-1 

i(4) 

i 

■ 

*1-1 

- a^*^ 
®1 

*1-1 

(7) 

‘i 

- 

*1-1 

• 

*1-1 

- a^’^ 
"l 

*1-1 

(7) 

,<2) 

A6) 

,(5) 

,(8) 

1 . 

1 

*i-1  • 

■ ®i- 

*1-1 

‘ «1 

*1-1 

(7) 

,(3) 

,(6) 

,(9) 

,(9) 

1 . 

1 

*i-1  • 

• *1 

*1-1 

■ «1 

*1-1 

.no  «'•’  - b'"’ 


(1^  ~ a 


for  1 < « <9 


- a 


(1) 

(1) 

^(2) 

(2) 

.(3) 

(3) 

'1 

^1-1  ■ 

«1 

^1-1  ■ 

■ »1  • 

M-1 

(4) 

(1) 

.(5) 

(2) 

(6) 

(3) 

1 

^1-1  ■ 

• *1 

^1-1  ■ 

■ «1  ^1-1 

C7> 

(1) 

(8) 

(2) 

C9) 

(3) 

1 

1^1-1  *1 

Yi-1  " Yi-1 

Kith  » d{"^  for  1 < ■ < 3 


Tht  arc  used  to  coaiputa  several  acre  functions  and  need  not  be 


stored  after  the  coaputatlon  of 
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.(1)  _ «C5>-(9)  . a(6)»(8) 

9^  * P^  P^  P^ 


.(6). (7)  . .<A)-(9) 


.(4)-(8)  , B<5)g(7) 
6.  . B.  B. 


.(4)  . .(3). (8)  ^ »(2)-C9) 

.(5)  . »<1)-(9)  _ -(3). (7) 


(2). (6)  .(3). (5) 


b:  b^' 


Br^V 


e<«>  = bP>b{^>  - Br>B!^> 

9,^’>  = 


u.  = ♦ e«>BP' 


The  and  are  used  to  compute  several  more  functions  and 

need  not  be  stored  after  the  computation  of 


,<1>  - ♦ .C4)^(4)  . .(7)^(7)., 

* (e^  ♦ 9^  ♦ 0^  c^  )/p^ 


<2)  . * .<7),(8)w 

* <9^  ♦ 9^  l/p^ 


The  values  of  X 


Must  be  stored,  as  they  are  used  in  the 


back  solution.  This  is 


